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ABSTRACT 

Context. This is the first paper of a series describing our measurement of weak tensing by large-scale structure, also termed "cosmic shear", 
using archival observations from the Advanced Camera for Surveys (ACS) on board the Hubble Space Telescope (HST). 
Aims. In this work we present results from a pilot study testing the capabilities of the ACS for cosmic shear measurements with early parallel 
observations and presenting a re-analysis of HST/ACS data from the GEMS survey and the GOODS observations of the Chandra Deep Field 
South (CDFS). 

Methods. We describe the data reduction and, in particular, a new correction scheme for the time-dependent ACS point- spread-function (PSF) 
based on observations of stellar fields. This is currently the only technique which takes the full time variation of the PSF between individual 
ACS exposures into account. We estimate that our PSF correction scheme reduces the systematic contribution to the shear correlation functions 
due to PSF distortions to < 2 x 10"^ for galaxy fields containing at least 10 stars, which corresponds to < 5% of the cosmological signal 
expected on scales of a single ACS field. 

Results. We perform a number of diagnostic tests indicating that the remaining level of systematics is consistent with zero for the GEMS 
and GOODS data confirming the success of our PSF correction scheme. For the parallel data we detect a low level of remaining systematics 
which we interpret to be caused by a lack of sufficient dithering of the data. Combining the shear estimate of the GEMS and GOODS 
observations using 96 galaxies arcmin"" with the photometric redshift catalogue of the GOODS-MUSIC sample, we determine a local 
single field estimate for the mass power spectrum normalisation crg coFS = 0.52+[J |j(stat) ± 0.07(sys) (68% confidence assuming Gaussian 
cosmic variance) at a fixed matter density Q.^ = 0.3 for a ACDM cosmology marginalising over the uncertainty of the Hubble parameter and 
the redshift distribution. We interpret this exceptionally low estimate to be due to a local under-density of the foreground structures in the CDFS. 
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Cosmic shear, the weak gravitational lensing effect of the large- 
scale structure, provides a powerful tool to constrain the to- 
tal matter power spectrum without any assumptions on the re- 
lation between luminous and dark matter. Due to the weak- 
ness of the effect, it is challenging to measure, with the first 
detec tions only reported six years ago (Bacon et al.l 200C ; 



KaiseretalJ I2OOOI; IVan Waerbeke et alj boOOt IWittman et al 



l2000i) . Since then cosmic shear has developed into a flourishing 
field of cosmology yielding not only constraints on the mat- 
ter content and the normalisation of the power spectrum 
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cently also on the equation of state parameter w of dark energ 



dJarvis et al.l2006HHoekstra et al.l2006l;ISemboloni et al.l200^ 
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for reviews. 

Due to the weakness of cosmological gravitational shear, 
proper correction for systematics, first of all the image point- 
spread-function (PSF) is indispensable. Within the Shear 
TEsting Programme (STEP) a number of algorithms to mea- 
sure the shear from faint and PSF-distorted galaxy images are 
cuiTently tested using a blind analysis of image simulations 
aiming to im prove and quantify the accuracy of the diffe rent 
methods; see iHevmans et al.l(l2006bl) : lMassev et al.l(l2007ah for 
first results. 

The majority of the previous cosmic shear measure- 
ments have been made with ground-based wide-field imag- 
ing data, mainly probing the matter power spectrum on lin- 
ear to moderately non-linear angular scales from several de- 
grees down to several arcminutes. In order to probe the highly 
non-linear power spectrum at arcminute and sub-arcminute 
scales, a high number density of usable background galax- 
ies is required. While deep ground-based surveys are typi- 
cally limited to ^ 30 galaxies/arcmin^ due to seeing, sig- 
nificantly higher number densities of resolved galaxies can 
be obtained from space-based images. Cosmic shear studies 
have already been ca r ried out with th e HST c ameras WFPC2 



dRhodes et al 



I2OOI'; Refregier et al. 2002; Casertano et al. 

2003h and STIS (IHammerle et al.. .2002; Rhodes et al.. .2004 ; 
Miralles et al. 20051) . With the installation of the Advanced 



Camera for Surveys (ACS) Wide-Field Channel (WFC) de- 
tector, a camera combining improved sensitivity (48% total 
throughput at 660 nm) and a relatively large fie ld-of-view 



(~ 3' 3 X 3'3) with good sa mpling (0'.'05 per pixel) (IFord et al 



2003t ISirianni et al.l 12004ft . the possibilities to measure weak 
lensing with HST have been further improved substantially. 
ACS has alread y been u sed for weak lensing measureme nts of 
galaxy clusters ( Jee et al...2005a.b..2006; Lombardi et al..i2005 



Clowe et al.ll2 006 HBradac et al.l2006tlLeonard et al.ll2007h and 



galaxy-galaxy lensing (IHevmans et al.l I2006at iGavazzi et al 
12007) . The first cosmic s h ear an a lysis with ACS was pre- 
sented by Hevmansetal.! (2005 ). H05 henceforth, for the 
GEMS survev dRix et all 12004*) a ~ 28' x 28' mosaic incor- 
porating the HST/ ACS GOODS observations of th e Chandra 
Deep Field Sott t h (CD FS, Gi avalisco et all 12004 . Recently 
Leauthaud et al.l (120071) and iMassev et al. I (12007b') presented a 



cosmological weak lensing analysis for the ACS COSMO^l 
field. 

In this work we present results from a pilot cosmic shear 
study using early data from the ACS Parallel Cosmic Shear 
Survey (proposals 9480, 9984; PI J. Rhodes). Parallel obser- 
vations provide many independent lines of sight reducing the 



impact of cosmic variance. With a separation of several ar- 
cminutes from the primary target (e.g. ~ 6' for WFPC2) they 
provide nearly random pointings for most classes of primary 
targets. This is important when measuring a statistical quantity 
like the shear. Still, for primary observations pointing at partic- 
ularly over-dense regions, such as galaxy clusters, a significant 
selection bias might be introduced as they influence the shear 
field many arcminutes around them, which has to be checked 
carefully. 

Analysing stellar fields we detect short-term variations of 
the ACS PSF, which are interpreted as focus changes due 
to thermal b r eathin g of the telescope ( see al so KristI 2003 



Rhodes et al.1 120051; [Anderson & Kind 120061; [Rhodes et al 



2007). Whereas earlier studies with other HST cameras as- 
sumed temporal stability of the PSF, a fully time-dependent 
PSF correction is required for ACS due to these detected varia- 
tions. Given that only a low number of stars are present in high 
galactic latitude ACS fields (~ 10 - 30), the correction cannot 
be determined from a simple interpolation across the field-of- 
view, but requires prior knowledge about possible PSF patterns. 
In this work we apply such a coiTection scheme, which is based 
on PSF models derived from stellar fields. Our method takes 
the full PSF variation between individual exposures into ac- 
count a nd can be applied for arbi trary dither patterns and ro- 
tations. Rhodes et al.l ( 2005 , 2007b propose a different correc- 
tion scheme, in which they fit co-added frames with theoretical 
single-focus PSF models created with a modified version of 

TinyTiiti. 
p — I 

H05 use a semi-time-dependent model to correct for the 
image PSF, based on the combined stars for each of the two 
GEMS observation epochs. The tests for systematics presented 
by H05 indicate zero contamination with systematics for the 
GEMS only data, but remaining systematics at small scales 
if the GOODS data are included. In order to test our fully 
time-dependent PSF correction we present a re-analysis of the 
GEMS and GOODS data, yielding a level of systematics con- 
sistent with zero for the combined dataset. We present a cos- 
mological parameter estimation from the re-analysed GEMS 
and GOODS data in this work, whereas a parameter estimation 
from the parallel data will be provided in a future paper on the 
basis of the complete ACS Parallel Cosmic Shear Survey. 

The paper is organised as follows; After summarising the 
weak lensing formalism applied in Sect.|2] we describe the data 
and data reduction in Sect.|3] We present our analysis of the 
ACS PSF and the correction scheme in Sect.|4] Next we elab- 
orate on the galaxy selection and determined redshift distri- 
bution (Sect.|5]l and compute several estimators for the shear 
and systematics in Sect.|6] After presenting the results of the 
cosmological re-analysis of the GEMS and GOODS data in 
Sect.|7] we conclude in Sect.|8] 

2. Method 

Cosmic shear provides a powerful tool to investigate the 3- 
dimensional power spectrum of matter fluctuations Ps through 
the observable shear field y - y\ + iyi induced by the tidal 



' |http : //www ■ astro . caltech . edu/''cosmos| 



^ |http : //www ■ stsci . edu/ sof tware/tinytim/tinytim. html 
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gravitational field (see iBartelmann & Schneideii uOOll for a 
broader introduction). They are related via the projected 2- 
dimensional shear (convergence) power spectrum 



4c4 



Jo 



dw- 



(1) 



where Ho is the Hubble parameter, the matter density pa- 
rameter, a the scale factor, w the comoving radial distance, Wh 
the comoving distance to the horizon, t the modulus of the 
wave vector, and /a-(w) denotes the comoving angular diame- 
ter distance. The source redshift distribution p„, determines the 
weighted lens efficiency factor 



g{w) = I dw /?„.(w )- 

Jw 



(2) 



see. e.g. lKaiseii (ll998l) : lSchneider et all (119981) 



2.1. Cosmic shear estimators 



Note, however, that the integral in (|9]l formally extends to in- 
finity. Thus, due to finite field size, real data require the sub- 
stitution of the measured ^-(ff) with theoretical predictions for 
large 9. 

This problem does not occur for the aperture mass statistics 



(ISchneider.1996.) 



(10) 



which is defined using the tangential and cross-components of 
the shear relative to the aperture centre ^ 



rt(0';^) = -^[r(eV^] 
rx(fl';^) = -3[r(0')e-2'^] 



(11) 

(12) 



with 0' - ^ - W - ^|(cos0 + isin0). For the axially- 
symmetric weight fun ction Qi'ff) we use a form proposed in 
Schneider etalTd 19981) 



In this work we measure the shear two-point correlation func- 
tions 



<rtrt)(0) = 



(rxTxXfl) = 



2f,- w,w;rx,i(-«^i) ■ jxjixj) 



(4) 



from galaxy pairs separated by = |jc, - Xj\, where the tangen- 
tial component ji and the 45 degree rotated cross-component 
7x of the shear relative to the separation vector are estimated 
from galaxy ellipticities, and w, denotes the weight of the rth 
galaxy. It is useful to consider the combinations 

^±ie)^{y,y,m±{y^y^m, (5) 

which are directly related to the convergence power spectrum 



27r Jo 



27r Jo 



(6) 



(7) 



w here J„ deno tes the n*-order Bessel function of the first kind. 



Crittenden et al.l (l2002h show that ^± can be decomposed into a 
curl-free 'E'-mode component ^^(0) and a curl 'B'-mode com- 
ponent ^^(0) as 



(8) 



6 (§ 



Q{ff) = — - 



H(6 - 1?) , 



(13) 



wh ere H(x) denotes the He aviside step function. 



Crittenden et al.l (l2002h show that Map purely measures the 



E-mode signal, whereas Mj_ contributes to the B-mode only. 

The dispersion of the aperture mass {MI^){9) is related to 
the convergence power spectrum by 



(14) 



with Wa^irj) - 576 14(77) rj^'^, and can in principle be computed 
by placing apertures on the shear field, yet in practice is prefer- 
entiall y calculated from the s hear two-point correla tion func- 
tions dCrittenden et al.l l2002t ISchneider et al.l 120021) to avoid 
problems with masked regions 



«p>(0) 



{Ml)(0) 



- r 

2 Jo 
2 Jo 



d^§ 



(15) 



(16) 



with T± given in lSchneider et al.l ( l2002h . 



2.2. Tine KSB formalism 



Kaiser. Squires. & BroadhurstI (119951) . iLuppino & Kaiser 
(I1997I) and lHoekstraet'air(ll998h (KSB-h) developed a 
formalism to estimate the reduced gravitational shear field 



with 



—^_{§)-n0' -^^-i-^)- (9) 

As weak gravitational lensing only contributes to the curl- 
free E-modes, such a decomposition provides an important test 
for the remaining contamination of the data with systematics. 



(17) 



from the observed images of background galaxies correcting 
for the smearing and distortion of the image PSF. In this for- 
malism the object ellipticity parameter 



— e\ + 162 — 



611-622+21612 
611 + 622 



(18) 
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is defined in terms of second-order brightness moments 

= J d^eWr^(\e\)6iejIie), i,je{l,2}, (19) use the approximations 



a function of the filter scale following lHoekstra et alj (Il998l) . 



For the calculation of P^^^ in (I2TI 1 and its inversion in (|24] | we 



where W,-, is a circular Gaussian weight function with filter 
scale Tg. The total response of a galaxy ellipticity to the reduced 
shear g and PSF effects is given by 



(20) 



with the intrinsic source ellipticity e^, the "pre-seeing" shear 
polarisability 



pg 



a/S ' 



psm f /psmt\-l pshtl 
^Q-y [V )y6 ^Sfi J 



(21) 



and the shear and smear polarisability tensors P^^ and P*™, 
which are calculated from higher -order brightness moments as 
detailed in lHoekstra et alJ (119981) . The anisotropy kernel q*{ff) 
describes the anisotropic component of the PSF and has to 
be measured from stellar images (denoted with the asterisk 
throughout this paper), which are not affected by gravitational 
shear and have e" - 0: 

ql^iP'^af^e;- (22) 
We define the anisotropy corrected ellipticity 

c'=^.-n>;^ (23) 

and the fully corrected ellipticity as 



„iso _ /pi>^.-l ani 



(24) 



which is an unbiased estimator for the reduced gravitational 
shear {e'^°) = g, assuming a random orientation of the intrin- 
sic ellipticity e^ For the weak distortions measured in cosmic 
shear <K 1 , and hence 



■ r- 



(25) 



The KSBh- formalism relies on the assumption that 
the image PSF can be described as a convolution of an 
isotropic part with a small anisotropy ker nel. Thus, it is ill- 
defined for several realistic PSF types (IKaisen l2000t) . be- 
ing of particular concern for diffraction limited space-based 
PSFs. Th is shortcoming incited the development of alternative 



Fiii^s. 1 n is snortcommg mcitea tne aeveiopment or alternative 
metho ds ( Rhodes et al. 2000t Kaiseil 2000t Bernstein & Jarvis 



2002'; 'Refregier & Bacon' '2003'; Massev & RefregieiJ 120051 
Kuiiken 2006 ; Nakajima & Bernstein 2007). Nevertheless 
Hoekstra et al.l (fl998l) demonstrated the applicability of the for- 
malism for HST/WFPC2 images, if the filter scale r„ used to 
measure stellar shapes is matched to the filter scale used for 
galaxy images 



In this pilot study we restrain the analysis to the Erben et al 



(1200 ih implementation of the most commonly used KSB-H for- 
malism. There are currently several independent KSB imple- 
mentations in use, which differ in the details of the compu- 
tation, yielding slightly different results (see Heymans et alJ 
2006bl for a comparison of several implementations). In partic- 
ular, our implementation interpolates between pixel positions 
for the calculation of Qij, P^^, and and measures all stellar 
quantities needed for the correction of the galaxy ellipticities as 



r/psm*\-l psh*l 
[V )yS ^dfi J 



Tr [P^h.j 
Tr [P™*] ' 



Tr [Ps] 



SaB , (26) 



as the trace-free pa rt of the tensor is much smaller than the trace 
dErben et al. 2001). To simplify the notation in the following 
sections we define 



Tr [P'^*] 
Trip™*] 



(27) 



We have tested this irnplementation with image simula- 
tions of the STEP projecjj. In the analysis of the first set 
of image si mulations (STEPl) we identified significant biases 
dHevmans et al..200 6b). which we eliminate with improved se- 
lection criteria (see Sect. l5.1b and the introduction of a shear 
calibration factor 



(28) 



with Ccai - 1/0.91. In a blind analysis of the second set 
of STEP image simulations (STEP2), which takes realistic 
ground-based PSFs and galaxy morphology into account, we 
find that the shear calibration of this improved method is on 
average accurate to the ~ 3% levefl. The method is capable 
to reduce the impact of the highly anisotropic ground-based 
PSFs which were analysed, to a level < 7 x lO""* dMassev et al. 
|2007a). We will also test this method on a third set of STEP 
simulations with realistic space-based PSFs (Rhodes et al. in 
prep.). Depending on the results we will judge whether the 
systematic accuracy will be sufficient for the complete ACS 
Parallel Survey or if a different technique will be required for 
the final cosmic shear analysis. For the GEMS Survey analysed 
in Sections|6]and|2l a ~ 3% calibration error is well within the 
statistical noise. In this work we use uniform weights w, - 1 
for all galaxies in order to keep the analysis as similar to our 
original STEP2 analysis as possible. 

3. Data 

The ACS Wide-Field-Channel detector (WFC) consists of two 
2k X 4k CCD chips with a pix el scale of ^05 y ielding a field- 
of-view (FOV) of ~ 3:3 X 33 dFord et aklbOOl . 

In this pilot project we use pure parallel ACS/WFC F775 W 
observations from HST proposal 9480 (PI J. Rhodes), denoted 
as the "parallel data" for the rest of this paper. 

For comparison we also apply our data reduction and anal- 
ysis to the combined F606W ACS/WFC observations of the 

GEMS field tRix et al 2004.) and the GOODS observations o f 

the Chandra Deep Field South (CDFS. lGiavaUsco et alj|2004 . 
A cosmic shear an alysis of this ~ 28' x 28' mosaic has already 
been presented by lH05l allowing us to compare the different 
correction schemes applied. 



' |http : //www ■ physics . ubc . ca/''heyman s/ ste p . html ] 
Note that detected dependencies of the shear estimate on size and 
magnitude will lead to slightly different uncertainties for different sur- 
veys. 
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Both datasets were taken within the first operational year 
of ACS (August 2002 to March 2003 for the parallel data; July 
2002 to February 2003 for the GEMS and GOODS observa- 
tions). Therefore these data enable us to test the feasibility of 
cosmic shear measurements with ACS at an early stage, when 
the charge-transfer-efRciency (CTE) has degraded only slight ly 
dRiess & Mackll2004l: lRiessll2004l: iMutchler & Siriannill2005h . 



3.1. The ACS parallel data 

The data analysed consist of 860 WFC exposures, which 
we associate to fields by joining exposures dithered by less 
than a quarter of the field-of-view observed with the guid- 
ing mode FINE_LOCK. In order to permit cosmic ray rejec- 
tion we only process associations containing at least two ex- 
posures. Furthermore, in this pilot study we only combine ex- 
posures observed within one visit and with the same role-angle 
in order to achieve maximal stability of the observing condi- 
tions. With these lim itations, which are similar to those used by 
Pirzkal et al.1 (12001 ) for the STIS Parallel Survey, we identify 
208 associated fields (including re-observations of the same 
field at different epochs), combining 835 exposures. 

For a weak lensing analysis, accurate guiding stability is 
desired in order to minimise variations of the PSF In case of 
parallel observations, differential velocity aberration between 
the primary and secondary instrument can lead to addit ional 
drifts during observations with the secondary instrument (ICox 
19971). In order to verify the guiding stability for each expo- 
sure we determine the size of the telemetry jitter-ball, which 
describes the deviation of the pointing from the nominal posi- 
tion. While the jitter-balls typically have shapes of moderately 
elliptical {{b/a} - 0.68) distributions with (FHWM) - 9.8 mas 
(0.196 WFC pixel), we have verified that FHWM < 20 mas 
(0.4 WFC pixel) and b/a > 0.4 for all selected exposures. 
Therefore the tracking accuracy is sufficiently good and ex- 
pected to affect the image PSF only slightly. Any residual im- 
pact on the PSF will be compensated by our PSF correction 
scheme, which explicitly allows for an additional ellipticity 
contribution due to jitter (see Sect. l4.4l l. 



3.1 .1 . Data reduction 

The data retrieve d was bias and flat-fiel d corrected. We use 
MultiDrizzl^ (Koekemoer et al. 2002) for the rejection of 
cosmic rays, the correction for geo metric camera dis tortions 
and differential velocity aberration ( Cox & Gillilandl 12002). 
and the co-addition of the exposures of one association. 

We refine relative shifts and rotations of the exposures 
by applying t he IRAF task geomap to matched windowed 
SExtractor (IBertin & ArnoutsI Il99 6) positions of compact 
sources detected in separately drizzled frames. For the star and 
galaxy fields selected for the analysis (see Sect. I3.1.2l i the me- 
dian shift refinement relative to the first exposure of an associ- 
ation is 0.17 WFC pixels, with 7.3% of the exposures requir- 
ing shifts larger than 0.5 WFC pixels. Refinements of rotations 
were in most cases negligible with a median of 1.6 x lO"'* de- 



grees corresponding to a displacement of ^ 0.008 WFC pixels 
near the edges of the FOV. Only in 1 .5% of the exposures rota- 
tion refinements exceeded 3 x 10"^ degrees, corresponding to 
displacements of ^ 0.15 WFC pixels. 

Deviat ing from the de fault parameters, we use the minmed 
algorithm (jPavlovsky et al..2006i) during the creation of the me- 
dian image as it is more efficient to reject cosmic rays for a 
low number of co-added exposures. For the cosmic ray masks 
we let rejected regions grow by one pixel in each direction 
(driz_cr_grow=3) in order to improve the rejection of neigh- 
bouring pixels affected due to charge diffusion and pixels with 
cosmic ray co-incidences in different exposures. 

Furthermore we use a finer pixel scale of 0'.'03 per pixel in 
combination with the SQUARE kernel for the final drizzle proce- 
dure in order to increase the resolution in the co-added image 
and reduce the impact of aliasing. For the default pixel scale 
(0'.'05 per pixel) resampling adds a strong artifical noise com- 
ponent to the shapes of un- and poorly resolved objects (alias- 
ing), which depends on the position of the object centre rela- 
tive to the pixel grid and most strongly affects the ei -ellipticity 
component. According to our testing with stellar field images, 
the GAUSSIAN kernel leads to even lower shape noise caused 
by aliasing. However, as it leads to stronger noise correlations 
between neighbouring pixels, we decided to use the SQUARE 
kernel for the analysis. 

Aliasing most strongly affects unresolved stars, which is 
critical if one aims to derive PSF models from a low number of 
stars in drizzled frames. However, since we determine our PSF 
model from undrizzled images (see Sect. 14. 4b . this does not af- 



fect our analysis. Consistent with the results from [Rhodes et al.l 
(2005, 2007) we find that a further reduction of the pixel scale 
does not further reduce the impact of aliasing significantly, 
while unnecessarily increasing the image file size. 

In this paper the term pixel refers to the scale of the drizzled 
images (0'.'03 per pixel) unless we explicitly allude to WFC 
pixels. 

3.1 .2. Field selection 

The 208 associations were all visually inspected. We discard in 
total 31 fields for the following reasons: 

- Fields which show a strong variation of the background in 
the pre-processed exposures (4 fields). 

- Fields containing galactic nebulae (10 fields). 

- Fields of significantly poorer image quality (6 fields). 

- Fields which contain a high number of saturated stars with 
extended diffraction spikes (5 fields). 

- Fields in M31 and M33 with a very high number density 
of stars, resulting in a strong crowding of the field, which 
makes them even unsuitable for star fields (4 fields). 

- Almost empty galactic fields affected by strong extinction 
(2 fields). 

Examples of the discarded fields are shown in Fig.[T] 

After this pre-selection, fields fulfilling the following crite- 



MultiPrizzle 2.7.0, |http : //stsdas . stsci . edu/multidrizzle/ j ria are selected as galaxy fields for the cosmic shear analysis: 
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Fig. 1. Examples for fields which were rejected by visual inspection for the following reasons (panels left to right): erroneous 
calibration, galactic nebula, many saturated stars, almost empty field. 



Table 1. Observation dates and position angles (ORIENTAT) of 
the ACSAVFC F606W GOODS/CDFS observations. 



Epoch 


Observation dates 


Position angle 


1 


2002-07-31-2002-08-04 


-112° 


2 


2002-09-19-2002-09-22 


-67° 


3 


2002-10-31-2002-11-03 


-22° 


4 


2002-12-19-2002-12-22 


23° 


5 


2003-02-01-2003-02-05 


68° 



- Fields have to be located at galactic latitudes \b\ > 25° in 
order to be affected only weakly by galactic extinction. 

- Only fields co-added from at least three individual expo- 
sures are used, facilitating sufficiently good cosmic ray re- 
jection. 

- Fields are required not to be dominated by a single object 
or stars resolved in a local group galaxy. 

- In the case of re-observations of the same field at differ- 
ent visits, the observation with the longest exposure time is 
used. 

55 independent fields fulfil these selection criteria. 
Additionally four fields with 20° < \b\ < 25° are included, 
which contain a high number density of galaxies indicating 
rather low extinction, making a total of 59 galaxy fields. This 
corresponds to 28.4% of the fields and 36.2% of the co-added 
exposures. 

All fields passing the preselection and containing at least 
300 stars are used as star fields for the PSF analysis (see 
Sect.m. These 61 fields consisting of 205 exposures amount 
to 29.3% of the fields and 24.5% of the co-added exposures. 

3.2. The GEMS+GOODS data 

The GEMS F606W data consist of 63 ACS/WFC tiles im- 
aged with three exposures of 720 to 762 seconds each. They 
are arranged around the ACS GOODS/CDFS observations, 
which have been imaged in five epochs with different posi- 
tion angles (see Tab.lTJ consisting of two exposures per tile and 
epoch with 480 to 520 seconds per exposure. In total the ACS 
GOODS/CDFS fieldis covered with 15 tiles during epochs 1, 3, 
and 5, whereas 16 tiles were used for epochs 2 and 4. H05 limit 
their analysis to the epoch 1 data. In order to reach a similar 



depth for the used GOODS and GEMS data we decided to com- 
bine the data of epoch 1 with either epoch 3 or 5 as they have an 
optimal overlap. The combination of epoch 1 and epoch 5 ex- 
posures is unproblematic. In contrast we find significant, FOV 
dependent residual shifts between matched object positions in 
exposures from epochs 1 and 3 after applying refined image 
shifts and rotations (Fig.|2]i. Possible interpretations for these 
remaining shifts are slight medium-term temporal changes in 
the ACS geometric distortion or a slightly imperfect treatment 
of the distortion correction in the MultiDrizzle version used. 



Pirzkal et al.l (120051) report similar effects for two epochs of 
Hubble Ultra Deep Field (UDF) data. As remaining shifts also 
occur for UDF images observed with position angles that are 
~ 90° apart, the MultiDrizzle interpretation might be more 
plausible. The largest residual shifts have a comparable mag- 
nitude of ~ 0.5 pixels for both the GOODS and UDF data. 
A combination of exposures with remaining shifts would in 
any case degrade the PSF of the combined image. Additionally 
central pixels of some stellar images could falsely be flagged 
as cosmic rays by MultiDrizzle. Therefore we only use the 
combined epoch 1 and 5 exposures for the cosmic shear analy- 
sis that follows. 

In order to investigate the ACS F606W PSF, we addition- 
ally analysed 184 archival F606W exposures of dense stellar 
fields containing at least 300 stars, which were observed be- 
tween July 2002 and July 2003. 

3.3. Catalogue creation 

We use SExtractor ( Bertin & Arnoutsll996 ) for the detection 
of objects and the Erben et al. (12001 ) implementation of the 
KSB formalism for shape measurements. We analyse the im- 
ages of galaxies in the combined drizzled images. However, for 
the time-dependent PSF correction described in Sect. 14.41 we 
additionally perform stellar shape measurements in the undriz- 
zled but cosmic ray-cleaned CO/? -images, which are also cre- 
ated by MultiDrizzle, and the drizzled uncombined frames 
(DRZ-images). 

The SExtractor object detection and deblending param- 
eters are summarised in Table|2l We use a rather low detec- 
tion threshold for the galaxies in order to minimise the im pact 
of PS F-based selection bias (Kaiser 2000; Bernstein & JarvisI 
l2002h . Spurious detections are later rejected with cuts in the 
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Fig. 2. Residual shifts [pixel] computed 
from windowed SExtractor positions of 
compact sources between epochs 3 and 1 
(left) and between epoch 5 and 1 (right) of 
the F606W ACS GOODS/CDFS observa- 
tions. For these plots compact objects from 
all 15 tiles are used, and residual shifts are 
averaged in bins of 700" pixels. For each tile 
the exposures of each epoch were drizzled 
onto one output pixel grid, with a common 
WCS per tile defined by epoch 1. Possible 
interpretations for the residual shifts in the 
left panel are slight temporal changes in the 
ACS geometric distortion or a slightly im- 
perfect treatment of the distortion correction 
in the MultiDrizzle version used. 



Table 2. Relevant parameters for the object detection with 
SExtractor for the galaxy fields and the star fields. Note that 
the number of pixels for a detection DETECTJIINAREA corre- 
sponds to the subpixel of the drizzled images except for the 
values in brackets, which are used for object detection in the 
undrizzled CO^-images. 



Parameter 


Galaxy fields 


Star fields 


BACK_TYPE 


AUTO 


MANUAL 


BACK_SIZE 


100 




BACK_FILTERSIZE 


5 




BACK_VALUE 




0.0 


DETECT_MINAREA 


16 


16(5) 


DETECT.THRESH 


1.5 


3(4) 


DEBLEND.NTHRESH 


16 


32 


DEBLEND.MINCONT 


0.05 


0.1 


FILTER_NAME 


gauss_2.5_5x5 (gauss_2.0_3x3) 



signal-to-noise ratio. We find that the deblending parameters 
applied perform well except for the case of spiral galaxies ex- 
tended by several arcseconds, for which substructure compo- 
nents are in some cases detected as separate objects. Thus, we 
mask these galaxies manually. If more than one object is de- 
tected within 1'.'2, only the brighter component is kept. We fur- 
thermore reject galaxies containing pixels with low values in 
the MultiDrizzle weight image (wmin = 100 sjl within their 
SExtractor isophotal area and also semi-automatically create 
masks to reject bright stars with diffraction spikes and extended 
image artifacts like ghost-images. 

We use different detection parameters for the star fields 
(see Table|2]i. Due to the increased detection threshold 
DETECT_THRESH, the object detection becomes less sensitive 
to the faint and extended stellar diffraction spikes, reducing the 
time needed for masking. 

We use the SExtractor FLUX_RADIUS parameter as 
Gaussian filter scale for the shape measurements of the 
galaxies. Here the integration is carried out to a radius of 3rg 
from the centroid. This truncation was introduced to speed up 



* The weight image pixel value corresponds to the effective expo- 
sure time contributing for the pixel, scaled with the relative area of 
output and input pixels. 



the algorithm and is justified due to the strong down-weighting 
of the outer regions in KSB. We also verified from the data 
that it does not bias the shape measurement. For the stellar 
D/?Z-images we repeat the shape measurements for 18 differ- 
ent filter scales ranging from 2.0 to 15 pixels, which are later 
matched to the filter scales of the galaxies. For larger filter 
scales we find that it is essential to continue the integration 
out to sufficiently large radii due to the wide diffraction wings 
of the PSF. Therefore, we employ a stellar integration limit of 
4.5 X FLUX_RADIUS* ^ 9 pixels. For the stellar shape measure- 
ments in the C(9/?-images we use a fixed Gaussian filter scale 
Tc, = 1.5 WFC pixels, which according to our testing roughly 
maximises the signal-to-noise of the stellar ellipticity measure- 
ment for most of the occuiTing PSF anisotropy patterns (see 
Sect.©. 

F or object se l ection we use the signal-to-noise definition 
from lErben etall (l200lb 



S/N = 



CTsky 



(29) 



which is based on the same filter function as the one used for 
shape measurements. In the computation of S/N we do not take 
the correlation of noise in adjacent pixels into account, which 
is created by drizzling. However, a correction for the noise in 
a large are a (e.g. the extent of a ga laxy) can be assessed from 
Eq. A19 in Casertano et al. ( 2000l) . which for our drizzle pa- 
rameters yields 



25/^ 
12 \ 9m 



(30) 



where 1x2 - cTsjiy is the single pixel background dispersion, 
while cr\ denotes the dispersion computed from areas of size 
n-p- (drizzled) pixels. The expression in squared brackets gives 
the correction factor to the area scaling expected for uncorre- 
lated noise. Using the effective area of the Gaussian weight 
function A - Inr-^ and m - yfA we estimate a noise correc- 
tion factor which increases from 1.86 for unresolved sources 
{rg ^2.1 pixels) to 2.05 for the largest galaxies considered. 
Hence, the true S/N will be lower than the directly computed 
value by this factor The cuts applied to the data refer to the 
directly computed value. 
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4. PSF analysis and correction 

Due to the low number of stars (~ 10 - 30) present in galaxy 
fields at high galactic latitudes we examined the ACS PSF from 
stellar fields (see Section I3.1.2l i containing ~ 300 - 20000 
stars. We do this analysis on the basis of single exposures in- 
stead of combined images, in order to optimally investigate 
possible temporal PSF variations. We investigate the PSF both 
in the undrizzled, but cosmic ray cleansed COR-images cre- 
ated by MultiDrizzle, and also the drizzled and cosmic ray 
cleansed single exposures (DRZ). Here we limit the discussion 
to the F775W data. Our analysis of the F606W PSF was per- 
formed in an identical fashion with only minor differences in 
the resulting PSF models. A det ailed KSB-H analysis of the 
F606W PSF can be found inlHOS . 



4.1. Star selection 

In the DRZ-images {COR-images) we select stars with 0.6 
pixel (0.45 WFC pixel) wide cuts in half-light radius rh 



(lErben et al.l 1200 Ih and cuts in the signal-to-noise ratio 



S/N > 40 (S/N > 30). We furthermore reject stars with satu- 
rated pixels using magnitude cuts, and, in the case of crowed 
fields, stars with a neighbour closer than 20 (10) pixels, which 
would otherwise affect the shape measurements for large r^. 

4.2. PSF anisotropy variation 

Investigating stellar fields we find that the stellar ellipticity 
e* and anisotropy kernel ^* vary smoothly across each WFC 
chip and can well be fit with third-order polynomials. Fig. [3] 
shows the FOV variation of e* for a 400 second stellar field 
exposure both for the undrizzled COR-image (left panel) and 
the drizzled and thus distortion corrected DRZ-image, where 
the middle panel corresponds to the PSF core measured with 
rg =2.4 pixels, whereas the right panel shows the PSF wings 
(rg = 10.0 pixels). The observed differences between the PSF 
core and wings, which mainly constitute in a stronger ellip- 
ticity for larger r^, underline the importance to measure stellar 



quantities a s a function of filter scale rg (see also lHoekstra et al 
19981: IhO?) 



In Fig.lH we compare the stellar ellipticity distribution in 
the COR-image and the DRZ-image for similar Gaussian filter 
scales of rg = 1.5 WFC pixels and rg - 2.4 pixels, both un- 
corrected and after the subtraction of a third-order polynomial 
model for each chip. Here drizzling with the SQUARE kernel in- 
creases the corrected ellipticity dispersion cr(ej"'*) by ^ 24% 
and thus decreases the accuracy of the ellipticity estimate. For 
the galaxy fields we therefore determine the PSF model from 
the undrizzled COR-images (see Sect. l4.4l i. Note that the stel- 
lar ellipticities in the COR-images (left panels in Fig.[3]and|4ll 
are created by the combined image PSF and geometric cam- 
era distortion, whereas the DRZ-image ellipticities correspond 
to pure image PSF. However, since the resulting pattern can 
in both cases well be fit with third-order polynomials, the cor- 
rected ellipticity dispersions are directly comparable. 

Note that we always plot the FOV variation in terms of 
e* in order to simplify the comparison to other publications. 
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Fig. 7. Histogram of the number of galaxy fields with A^star 
selected stars in the co-added images for the Parallel Survey 
(dashed fine) and the GEMS-hGOODS data (solid line). 



However, for the actual correction scheme we employ fits of q*^ 
defined in (l22T l due to a slight PSF width variation leading to a 
variation of f™* (see Sect.|43]l. 

Comparing stellar field exposures observed at different 
epochs, we detect significant temporal variations of the PSF 
anisotropy already within on e orbit . Time v ariations of the AC S 
PSF were also reported by 'Krist' ('2003'); 'jee et al.' ('2005a''b', 
1 2006); H05; Rhodes et al. (2005, 2007) and Anderson & King 
(120061) . and are expected to be caus ed by focus changes due to 



thermal breathing of the telescope. iKristI (12003) illustrates the 
variation of PSF ellipticity induced by astigmatism, which in- 
creases for larger focus offsets and changes orientation by 90° 
when passing from negative to positive offsets. This behaviour 
is approximately reproduced in Fig. |5] showing polynomial fits 
to stellar ellipticities in two series of subsequent exposures. 

4.3. PSF widtii variation 

Additional to the PSF ellipticity variation we also detect time 
and FOV variations of the PSF width. Fig.|6] shows the FOV 
dependence of the stellar half-light radius rh for three different 
exposures. Among all F775W stellar field exposures the aver- 
age half-light radius varied between 1.89 < rh < 2.07, with an 
average FOV variation cr(rh) - 0.085. We find that the vari- 
ation of the stellar quantity T* needed for the PSF correction 
of the galaxy ellipticities ( l27T i closely follows the variation of 
rh and can well be fitted with fifth-order polynomials in each 
chip. For a further discussion of the PSF width variation see 
lKrisd(l2003h . 



4.4. PSF correction scheme 

In order to correct for the detected temporal PSF variations us- 
ing the low number of stars present in most galaxy fields (see 
Fig.|7]l, we apply a new correction scheme, in which we de- 
termine the best-fitting stellar field PSF model for each galaxy 
field exposure separately. 
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Fig. 3. Stellar "whisker plots" for an example F775W stellar field exposure. Each whisker represents a stellar ellipticity. The 
left panel shows e* in the undrizzled COR-image measured with rg = 1.5 WFC pixels (PSF core). The middle and right panels 
correspond to the drizzled DRZ-image showing the PSF core (rg = 2.4 pixels, middle) and the PSF wings (rg - 10.0 pixels, 
right). The fit to the ellipticities in the middle panel is shown in the lower right panel of Fig.|5] 
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Fig. 4. Stellar ellipticity dis- 
tribution (PSF core) for an 
example F775W stellar field 
exposure measured in the 
undrizzled COiJ-image (left) 
and the drizzled DRZ-image 
(right). The open circles rep- 
resent the uncorrected el- 
lipticities e*, whereas the 
black points show the ellip- 
ticities e™'* corrected with 
a third-order polynomial for 
each chip. In the right panel 
cr{e'f"') is significantly in- 
creased, which is a result of 
the re-sampling in the drizzle 
algorithm. In these plots out- 
liers have been rejected at the 
3(T level. 



4.4.1. Description of the algoritlim 

Due to the low number of stars present in galaxy fields, we 
require a PSF fitting method with as few free parameters as 
possible, excluding the possibility to use for example a direct 
polynomial interpolation. As the main PSF determining factor 
is the focus position, we expect a nearly 1 -parameter family of 
PSF patterns. With the high number of stellar field exposures 
analysed A^sf = 205 for F775W and N,f = 184 for F606W, we 
have a nearly continuous database of the varying PSF patterns. 
This database consists of well constrained third- or fifth-order 
polynomial fits to qaix, y, rg) and T(x, y, r„) for numerous val- 
ues of rg, both for the COR- and D7?Z-images. In this section 
we omit the asterisk when we refer to these polynomial fits 
derived from the stellar fields in order to allow for a clear dis- 
tinction to ql, measured from the stars in the galaxy fields. 

Given the noisier e* and ^* measurement in drizzled im- 
ages (see Sect. 14. 21 ). we estimate the PSF correction for a galaxy 
field from the stellar images in each C(9/?-exposure of the 
galaxy field. However, we apply the corresponding DRZ-image 



PSF model as galaxy shapes are also measured on the drizzled 
co-added image. 

In order to determine the correction for a co-added galaxy 
field, we fit the measured ga'^'^^ir^ - 1.5) of the Astars.t stars 
present in galaxy field exposure k with the stellar field PSF 
models q'^'^^(x,y, r^), with j e 1, ...,Asf and identify the best 
fitting stellar field exposure jk with minimal 

- Z ^'•^ = 1-5) - ^"f (-"3'" 1-5)1 ■ (31) 

(=1 

Here we choose the Gaussian window scale rg = 1 .5 WFC pix- 
els to maximise the signal-to-noise in the shape measurement 
(see Sect. l3.3l l. In this fit we reject outliers at the 2.5cr level to 
ensure that stars in the galaxy field with noisy ellipticity esti- 
mates do not dominate the fit. 

Having found the "most similar" (best-fitting) COR-PSF 
model jit for each galaxy field exposure, we next have to match 
the coordinate systems of the corresponding Di?Z-image and 
the co-added galaxy field. This is necessary, as the single DRZ- 
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Fig. 5. Third-order polynomial fits to stellar ellipticities in the DRZ-images of two series of subsequent exposures. The 400 second 
exposures were taken on 2002-08-28 (upper panels) and 2002-08-17 (lower panels), where the time indicated corresponds to the 
middle of the exposure (UT). The variations are interpreted as thermal breathing of the telescope. The upper right and lower 
left plots are near the optimal focus position, whereas the other exposures represent positive focus offsets (upper left panel) or 
negative focus offsets (lower right panel). 
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Fig. 6. Field-of-view variation of the stellar half-light radius in the D/?Z-images of three F775W stellar field exposures with 
positive focus offset (left), near optimal focus (middle) and negative focus offset (right). For these plots rt was averaged within 
bins of 300 x 300 pixels. Bins without any stars show the average value. 



images used to create the PSF models are always drizzled with- 
out extra shifts in the default orientation of the camera, whereas 
the galaxy field exposures are aligned by MultiDrizzle ac- 
cording to their dither position. For this we trace the position 



of each object in the co-added galaxy field back to the position 
it would have in the single drizzled exposure k without shift 
and rotation. Let and (xo, yo)k denote the rotation and shift 
appUed by MultiDrizzle for exposure k. For a galaxy with 
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coordinates (x, y) in the co-added image we then compute the 
"single Di?Z"-coordinates 



cos sin 



■ sm( 



cos (f>k \y 



■yo,k 



and the PSF model 



l\<pt 



T^^ix,y,r,) 



(32) 



(33) 
(34) 



where we denote the components of as q^^- 

In order to estimate the combined PSF model for the co- 
added galaxy image, we then compute the exposure time tt- 
weighted average 



?°,comb(-^'3','-g) 



Y^hq^f(x,y,r^)Ak / ^ f^A, (35) 

k ) k 

Y^hT^^(x,y,r^)Au /^f^A, (36) 



of all shifted and rotated single exposure models, with Aj = 1 
if the galaxy is located within the chip boundaries for exposure 
k and A^^ = otherwise. 

Another factor which is expected to influence the image 
PSF besides focus changes are jitter variations created by track- 
ing inaccuracies (Sect.|3]l. To take those into account we fit an 
additional, position-independent jitter term ^°(rg). We already 
take this constant into account while fitting the galaxy field 
stars with the stellar field models to ensure that a large jitter 
term does not bias the identification of the best-fitting star field. 
Yet, as the number of stars with sufficient signal-to-noise is 
higher in the co-added image and since only the combined jitter 
effect averaged over all exposures is relevant for the analysis, 
we re-determine the jitter term in the co-added drizzled image 
after subtraction of the combined PSF model q^^^ij'.x,y,rg). 
The final PSF model used for the correction of the galaxies is 
then given by 



1^Z.M^y^ - ?°comb(^-3', rg) + ql(r,) 



(37) 



andrD^Z(x,y,rg). 

Note that this correction scheme assumes that the PSF 
model quantities q^^ix, y, rg) and T^^^(x, y, r^) determined 
for each galaxy field exposure can directly be averaged to deter- 
mine the correction for the co-added image. While only bright- 
ness moments add exactly linearly, this computation simplify- 
ing approach is still justified, as both the PSF size variation 
and the absolute value of the stellar ellipticities are small (see 
Sections [4.21 and l43T l. Computing the flux-normalised trace of 
the stellar second brightness moments 



Qu + Qi2 

FLUX* 



(38) 



for all stars in the F775W stellar field exposures with fixed 
Tg = 2.4 pixels, we find that Q has a relative variation of 3% 
only (Icr). Therefore we can well neglect non-linear terms in- 
duced by the denominator in (fTST l. The same holds for non- 
linear contributions of P^"* and T*, which show Icr- variations 



only by 6% (Tr[P™*]) and 5% (r*). As a (very good) 
first-order approximation we can therefore simply average 
q^^{x, y, Tg) and T^^^{x, y, r^) linearly, as also demonstrated 
in Fig.|8] where we compare the exposure time-weighted aver- 
age of the quantities measured from stars in individual drizzled 
frames to the value measured in the co-added image. 



4.4.2. Test with star fields 

In order to estimate the accuracy of our fitting scheme we test 
it on all co-added stellar field images. For each stellar field we 
randomly select subsets of A^star stars from the COR-images and 
the co-added image simulating the low number of stars present 
in galaxy fields. This subset of stars is used to derive the PSF 
model as described in Sect. l4.4.ri which we then apply to the 
entirety of stars in the co-added image. For the fitting of a par- 
ticular stellar field exposure, we ignore its own entry in the PSF 
model database and only consider the remaining models. The 
strength of any coherent pattern left in the stellar ellipticities 
after model subtraction provides an estimate of the method's 
accuracy. In order to determine the actual impact of the remain- 
ing PSF anisotropy on the cosmic shear estimate, one has to 
consider that although galaxy ellipticities are on the one hand 
less affected by PSF anisotropy than stars, they are additionally 
scaled with the correction (I23I24I I. We thus "transform" the 
remaining stellar ellipticity into a corrected galaxy ellipticity 
(see e.g. ,Hoekstraii2004) 



2Ccal 
gal 



gal 

TrP™*(rg) 



4('-g)-nAgai<omi(^>3'>'-g) 



,(39) 



where we randomly assign to each star the value of ff^j, 
and rg from one of the parallel data galaxies used for 
the cosmic shear analysis (see Sect.|6]l. Fig.|9] shows the esti- 
mate of the two-point correlation functions of e*;'^° averaged 
over all star fields and 30 randomisations for different num- 
bers of random stars A^star- This plot indicates that already for 
A^star = 10 stars present in a galaxy field the contribution of re- 
maining PSF anisotropy is expected to be reduced to a level 
<^r/x°^*/x°> < 2 X 10"^ corresponding to ^ 1 - 5% of the cos- 
mological shear correlation function expected on scales probed 
by a single ACS field. Since all of the examined parallel fields 
and the large majority of the GEMS-hGOODS fields contain 
more than 10 stars (see Fig.|7]l, we are confident that the sys- 
tematic accuracy of this fitting technique will be sufficient also 
for the complete ACS parallel data. 

The further reduction of the remaining systematic signal 
for larger A^star shows that the accuracy is mainly limited by the 
number of available stars and not by a too narrow coverage of 
our PSF database or the linear averaging of q^^^{x,y, r„). 

For comparison we also plot in Fig.|9]the correlation func- 
tions calculated from the f ^-scaled, but not anisotropy cor- 
rected stellar ellipticity, which for larger scales is of the same 
order of magnitude as the expected shear signal. Note that the 
plotted values depend on the selection criteria for the galax- 
ies (see Sect. 15. II ). Particularly, the inclusion of smaller, less 
resolved galaxies would increase both the corrected and uncor- 
rected signal. Additionally, it is assumed that the distribution 
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Fig. 8. Comparison of the stellar quantities ^* and T* measured on individual stars with rg — 2.4 pixels from a co-added stellar 
field to the same quantities computed as an exposure time-weighted average of the estimates in the single D/?Z-images. The good 
and unbiased agreement justifies the direct use of these quantities in the PSF correction scheme without the need to work on 
individual moments. The plotted points coiTespond to the three stellar field exposures shown in the bottom of Fig.|5] Note the 
larger scatter for q[ compared to q'2 which is mainly due to the noise created by re-sampling. 



of PSFs occurring is the same for the star and galaxy fields. For 
more homogeneous surveys (e.g. the GEMS+GOODS data) 
one might expect that more similar PSFs occur more frequently 
than for the quasi random parallel star fields, for which the 
stellar correlation function is expected to partially cancel out. 
Thus, we also plot the one sigma upper limit of the unconnected 
correlation function in Fig.|9l which might be a more realistic 
estimate for the uncorrected systematic signal for such surveys. 



4.4.3. Discussion of the algorithm 

The applicability of the proposed algorithm relies on the as- 
sumption that the stellar fields densely cover the parameter 
space of PSF patterns occurring in the galaxy fields. This is 
likely to be the case if 



1 . both datasets roughly cover the same time span, 

2. the number of star field exposures is sufficiently large, 

3. and no significant additional random component occurs be- 
sides the constant jitter offset that we have considered. 



For both the F606W and F775W data (1.) is fulfilled and from 
the ensemble of observed stellar field PSFs we are confident 
that (2.) and (3.) are also well satisfied. This is also confirmed 
by the test presented in Sect. 14.4.21 Yet, the reader should note 
that datasets might exist for which conditions (1.) to (3.) are 
not well fulfilled, e.g. due to observations in a rarely used fil- 
ter with only a low number of observed stellar fields. In such a 
case the described algor ithm might be adjus ted using a princi- 
pal component analysis ("jarvis & Jainl2004 ) or theoretical PSF 
models (Rhodes et al. 2005, 200_^ Note that the differences in 
the observed PSFs are interpreted to be mainly driven by differ- 
ent focus offsets. However, the suggested algorithm will work 
just as well if further factors play a role, as long as sufficient 
stellar field exposures are available. 



4.4.4. Advantages of our PSF correction scheme 

Finally we want to summarise ffie advantages our method pro- 
vides for the high demands of a cosmological weak lensing 
analysis on accurate PSF correction: 

1 . Our technique deals very well with the low number of stars 
present in typical high galactic latitude fields, which in- 
hibits direct interpolation across the field-of-view. 

2. The ACS PSF shows substantial variation already between 
consecutive exposures (see Fig.|5]l, which is adequately 
taken into account in our technique. When exposures from 
different focus position s are combine d , a sin gle-focus PSF 
model, as e.g. used by [Rhodes et alJ (l2007h . is no longer 
guaranteed to be a good description for the co-added im- 
age. 

3. Our PSF models are based on actually observed stellar 
fields and are thus not affected by possible limitations of 
a theoretical PSF model. 

4. We determine the PSF fits in the undrizzled C(9/?-images, 
which excludes any impact from additional shape noise in- 
troduced by re-sampling. 

5. The algorithm is applicable for arbitrary dither patterns and 
rotations, and can easily be adapted for offier weak lensing 
techniques (e.g. Nakajima et al. in prep.). 

5. Galaxy catalogue and redshift distribution 

5.1. Galaxy selection 

We select galaxies with cuts in the signal-to-noise ratio 
S/N > 4, half-light radius 2.8 < rh < 15 pixels, corrected 
galaxy ellipticity |e™| < 2.0, and TrP^/ 2 > 0.1. From the anal- 
ysis of the STEPl image simulations (iHevmans et al. I l2006bh 
we find no indications for a significant bias in the shear esti- 
mate introduced by these conservative cuts for |e'^°| and Trf ^. 
However, due to a detected coiTelation of the shear estimate 
both with rh and the SExtractor FLUX_RADIUS, cuts in these 
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Fig. 9. Estimate for the PSF fitting accuracy: In order to sim- 
ulate the low number of stars in galaxy fields, the PSF correc- 
tion technique was applied to the 61 parallel data star fields, 
from which only small random subsets of A^stai stars were 
used to determine the fit. We plot the correlation functions 
{e*/'°e*'''°} (left) and <e*'''°e*'''°> (right) of the "transformed" 
and corrected stellar ellipticity e*''^° (39[ . which accounts for 
the different susceptibility of stars and galaxies to PSF effects. 
The numbers A^star = (5,10, 20, 50) indicate the number of ran- 
dom stars used in each subset. Note that the uncorrected PSF 
signal computed from the transformed but not anisotropy cor- 
rected ellipticity (nocor) and its Icr upper limit (nocor-i-1 cr) 
are only slightly lower than ACDM predictions for the cosmo- 
logical lensing signal shown by the dashed-dotted curves for 
o-a =0.7, zm = 134. 



quantities introduce a significant selection bias. For the anal- 
ysis of the STEP2 image simulations we therefore chose r^ - 
cuts closely above the stellar sequence jMassev et al. 2007ah . 
Yet, as the magnitude-size relation is very different for ground- 
and space-based images, we expect that a cut at larger rh will 
introduce a smaller shear selection bias for space-based im- 
ages. In Fig. [To] we plot the rh -magnitude distribution of the 
objects in the F775W galaxy fields after a cut S/N > 4 was 
applied. Considering the PSF size variation (Sect. l4.3T l and in- 
creased noise in the rj, measurement for faint objects, stars can 
reliably be rejected with a cut > 2.4 pixels. With the cuts 
in le"*"!, and Trf* applied, increasing the size cut to > 2.8 
pixels rejects only 6.1% of the remaining galaxies. As these 
galaxies are most affected by the PSF, and considering the pos- 
sible limitations for the application of the KSB formalism for 
a diffraction limited PSF (Sect.|2]l, we decided to use the more 
resolved galaxies with > 2.8 pixels. We plan to investigate 
whether this introduces a significant shear selection bias with 
the space-based STEP3 simulations (Rhodes et al. in prep.). 

With these cuts we select in total 39898 (77749) galax- 
ies corresponding to an average galaxy number density 




6 8 
rh [pixels] 

Fig. 10. rh-magnitude distribution of the Parallel data F775W 
objects after applying a cut S/N > 4. The vertical lines indicate 
two different cuts for the galaxy selection: Although a cut rh > 
2.4 pixels is sufficient to reliably exclude stars, we additionally 
reject very small galaxies (2.4 pixels < rh < 2.8 pixels), which 
are most strongly affected by the PSF. 

of 63 arcmin"^ (96 arcmin"^) for the parallel F775W fields 
(GEMS+GOODS F606W tiles) with a corrected ellipticity dis- 
persion cr(ccaie™) = 0.32 (0.33) for each component. 

Ih05' found that the faintest galaxies in their catalogue were 
very noisy diluting the shear signal. Therefore they use a con- 
servative rejection of faint galaxies (meoe < 27.0, SNR > 15) 
leading to a lower galaxy number density of ~ 60 arcmin"^ for 
the GEMS and GOODS F606W data. For our primary anal- 
ysis we use a rather low cut S/N > 4 (see above) to be con- 
sistent with our analysis of the STEP simulations. In order to 
assess the impact of the faintest galaxies and ease the compari- 
son to H05, we repeat the cosmological parameter estimation in 
Sect.|7]with more conservative cuts S/N > 5, m606 < 27.0 lead- 
ing to a number density of 72 arcmin"^, which is roughly com- 
parable to the value found by H05i given the deeper combined 
GOODS images in our analysis. 

We plot the average galaxy number density as a function 
of exposure time for the different datasets in Fig.[TT] indicat- 
ing that F606W is more efficient than F775W in terms of the 
average galaxy number density. However one should keep in 
mind that the parallel fields are subject to varying extinction 
and more inhomogeneous data quality. 

For the GEMS-hGOODS tiles we rotate the galaxy elliptic- 
ities to a common coordinate system and reject double detec- 
tions in overlapping regions which leaves 71682 galaxies for 
S/N > 4 and 53447 galaxies for S/N > 5, mm < 27.0. 



5.2. Comparison of shear catalogues 

In regions where different GEMS and GOODS tiles overlap, we 
have two independent shear estimates from the same galaxies 
with different noise realisations corrected for different PSF pat- 
terns. This provides us with a good consistency check for our 
shear pipeline. We compare the two shear estimates in the left 
panel of Fig. [12] Although there is a large scatter created by the 
faint galaxies, which are strongly affected by noise, the shear 
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Fig. 12. Comparison of the shear 
estimates between overlapping 
ACS tiles (left) and between the 
|h05 and our catalogue (right). The 
grey-scale indicates the number of 
galaxies. Note the slight difference 
in the shear calibration between the 
two pipelines (~ 3.3%). In the left 
panel galaxies from different noise 
realisations are compared, leading 
to the larger scatter. The solid line 
shows a 1 : 1 relation. 



5.3. Redshift distribution 
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Fig. 11. Number density of selected galaxies (S/N > 4) for the 
parallel data F775W fields and the GEMS/GOODS F606W 
tiles as a function of exposure time. 



estimates agree very well on average confirming the reliability 
of the pipeline. 



Additionally, we match our shear catalogue to the |H05| cat- 
alogue, which stems from an independent data reduction and 
weak lensing pipeline, and compare the two shear estimates in 
the right panel of Fig. [12] Overall there is good agreement be- 
tween the two pipelines with a slight difference in the shear cal- 
ibration, where our shear estimate is in average larger by 3.3%. 
This is also consistent with results of the STEP project given 
a 3% under-estimation of the shear for the Heymans pipeline 
in STEPl (jneymans et al. 2()06b) and an error of the aver- 
age shear calibrati on consistent with ze ro for the Schrabback 
pipehne in STEP2 dMassev et al.ll2007ah . The slightly diff^erent 
results for the two KSB-H pipeUnes are likely to be caused by 
the shear calibration factor used in our pipeline and the differ- 
ent treatment of measuring shapes from pixelised data, where 
we interpolate across pixel s while H05 evaluate t he integrals at 
the pixel centres. See also iMassev et al.l (l2007al) for a discus- 
sion of the impact of pixelisation based on the STEP2 results. 
For the GEMS and GOODS data a shear calibration error of 
~ 3% is well within the statistical noise. 



In order to estimate cosmological parameters from cosmic 
shear data, accurate knowledge of the source redshift distri- 
bution is required. This is of particular concern if the red- 
shift distribution is constrained from external fields (see e.g. 
Ivan Waerbeke et al. ||200^ iHuterer et all l2006h . However, as 
the Chandra Deep Field South has been observed with sev- 
eral instruments including infrared observations, accurate pho- 
tometric redshifts can directly be obtained for a significant 
fraction of the galaxies without the need for external calibra- 
tion. In this work we use the photometric reds hift catalogue 



of the GOODS-MUSIC sample presented by iGrazian et al 



(120061) . This catalogue combines the F435W, F606W, F775W 



and F850LP ACS GOODS/CDFS images (iGiavalisco et al 



20041) . the JHKs VLT data (Vandame et al. in prep.), the 



Spitzer data provided by the IRAC instrument at 3.6, 4.5, 
5.8, and 8.0 fim (Dickinson et al. in prep.), and f/-band 
data from the MPG/ESO 2.2m and VLT-VIMOS. Additionally 
the GOODS-MUSIC catalogue c ontains spectrosc opic data 



Wolf et al.! 2001; Bunker et al 



from several surveys (ICristiani et a l. 2000; Croom et al 



2003; Dickinson et al 



Le Fevre et al.i|2004t jSzok olv et al .li2004i; iStanwav et al 



Strolger et al 



l2004t [ van 
a et al.l I2OO: 



der Wei et al. 
20051) 




Mignoli e t al.l 
which are also compiled in a 



20051; Vanzella et 

MasteJ!! catalogue by the ESO-GOODS team. We match the 
GOODS-MUSIC catalogue to our filtered galaxy shear cata- 
logue, yielding in total 8469 galaxies with a photometric red- 
shift estimate, including 408 galaxies with additional spectro- 
scopic redshifts and a redshift quality flag qz < 2. In the 
area covered by the GOODS-MUSIC catalogue 95.0% of the 
galaxies in our shear catalogue with mf,o6 < 26.25 have a 
redshift estimate, and only for fainter m agnitudes sub s tantia l 
redshift incompleteness occurs (Fig. [T3] l. lOrazian et al. (2006) 
estimate the photometric redshifts errors from the absolute 
scatter between photometric and spectroscopic redshifts to be 
<|Az/(l -Hz)|) = 0.045. 



' |http : //www ■ eso . org/ science/goods/spectroscopy/CDFS_Mastercat/| 
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Fig. 13. Number of selected GOODS-CDFS galaxies as a func- 
tion of m606- The solid line corresponds to galaxies for which 
spectroscopic or photomet ric redshift are available from the 
GOODS-MUSIC sample dOrazian et al] booel) . whereas the 
dotted line shows galaxies in the shear catalogue without red- 
shift estimate located in the same area. 



In cosmic shear studies the redshift distribution is often 
parametrised as 
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Fig. 14. Redshift distribution for the matched shear catalogue 
galaxies with redshift estimate from the GOODS-MUSIC sam- 
ple in the magnitude range 21.75 < meoe < 26.25 (solid line 
histogram). The dashed curve shows the reconstructed red- 
shift distribution N<p{z) for these galaxies using the best fitting 
values for (a, = (0.563, 1.716,0.299,0.310). The dotted 
curve was computed for fixed (a,yS) = (2, 1.5). 



Piz) 



/3 



— exp 



ZO 



(40) 



Brainerd et al 
2006') 



19961: ISemboloni et"ai] l2006t 



Hoekstra et al.l |2006'). In order to extrapolate the redshift 
distribution for the faint and redshift incomplete magnitudes 
we consider p(z) = piz,in(,o(,) and assume a linear relation 
between the magnitude meoe and the median redshift Zm of an 
ensemble of galaxies with this magnitude 

Zm = '-Z0 = fl('W606 - 22) + b, (41) 

where r(a,fi) is calculated from numerical integration of ( |40] l. 
For a single galaxy of magnitude meoe, (|40] | corresponds to 
the redshift probability distribution given the parameter set 
(a,/3, a, b). Thus, we can constrain these parameters via a max- 
imum likelihood analysis, for which we marginalise over the 
photometric redshift errors Az. The total redshift distribution 
of the survey with galaxies is then constructed as 



0(z) = 



N 



(42) 



Note that this approach is similar to the one used bv lH05L but 
does not require magnitude or redshift binning. 

For the maximum likelihood analysis we apply the CERN 
Program Library MINUItQ and use all galaxies with redshift es- 
timates in the magnitude range 21.75 < meoe < 26.25. Varying 
all four parameters {a,/3,a,b) we find the best fitting param- 
eter combination (a,/S,a,b) = (0.563, 1.716,0.299,0.310), for 
which Zm - 0.7477zo. In order to estimate the fit accuracy, we 
fix a and /3 to the best fitting values and identify the 68% (95%) 



Ihttp : //wwwasdoc . web . cern . ch/wwwasdoc/minuit/ I 



confidence intervals for a and b: a 



0.310 



+0.018(0.037) 



^^00+0.006(0.013) , 
^■^^^-0.007(0.014)' 



-0.017(0.033)' 

Using these parameter estimates, we reconstruct the 
redshift distribution of the galaxies used for the fitting 
(Fig.[T4b. The reconstruction fits the actual redshift dis- 
tribution very well except for a prominent galaxy over- 
density at z ^ 0.7 and an under-density at z > 1.5, 
which are known large-scale structure features of the field 
Gillietal. 2003, 2005; Szokolv et al. 2004; LeFevreetal 
20041; IXdami et all 120051; IVanzella et al] I2OO6I; iGrazian et al 



I). Yet, given that the reconstruction and the photomet- 
ric redshift distribution have almost identical average redshifts 
(2recon(fit sample)) = 1.39, (Zphoto(fit sample)) = 1.41, we esti- 
mate that the impact of the large-scale structure on the cosmic 
shear estimate via the source redshift distribution will be small. 
However, the large-scale structure significantly influences 
the estimate of the median redshift Zni,recon(fit sample) = 1.23, 
Zm,photo(fit sample) = 1.10. Thus, a redshift distribution deter- 
mined from the computed median redshift of the galaxies 
would most likely be biased to too low redshifts. Note that in 
Fig.[T4]the reconstruction falls off slower for high z than the ac- 
tual distribution of the data. To exclude a possible bias we thus 
always truncate the high redshift tail for z > 4.5. 

For comparison we also determine a reconstruction from 
the best fitting values for {a,b) with fixed values (a,/3) - 
(2, 1.5), w hich are sometimes used in the literature (e.g. 
Baugh & Ef stathiou 1994; H05). While they seem to pro- 



vide a good param etrisation for shallower surveys (see e.g. 



Brown et alj|2003h . they lead to a distribution that is too nar- 



rowly peaked with a maximum at too high redshifts for the deep 
GEMS and GOODS data (Fig.IH). 
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Fig. 15. Median redshift of the matched galaxies in mgoe 
bins computed directly from the data (thin crosses) and de- 
termined from a maximum likelihood fit for Zm with fixed 
ia,/3) = (0.563, 1.716) (triangles), with eiTors-bars indicating 
the eiTor of the mean or the Icr confidence region, respectively. 
The solid line corresponds to the best fitting parameters of the 
joint UkeUh ood fi t, whereas the dashed line shows the fit de- 
termined by 'HOS' for the magnitude range 21.8 < »J606 < 24.4. 
Note that a large-scale structure peak at z ^ 0.7 induces both 
the flatter slope for the directly computed Zm for meoe ^ 24.7 
and the increased spread for the fitted points for meoe 5 23.3. 
For mgoe Z 26.25 substantial redshift incompleteness occurs. 
For comparison we also plot the directly comput ed median 
photometric redshift from the HUDF (open circles, ICoe et al. 
[2006.) . 

A maximum likelihood analysis can only yield reasonable 
parameter constraints if the model is a good description of the 
data. To test our assumption of a linear behaviour in i4T[ . we 
bin the matched galaxies in redshift magnitude bins and de- 
termine a single Zm for each bin using an additional likeli- 
hood fit with fixed (a,/3) = (0.563, 1.716), see Fig.lB] A linear 
ZmCwJeoe) description is indeed in excellent agreement with the 
data in the magnitude range used for the joint fit. Only at the 
bright end the large-scale structure peak at z ^ 0.7 induces an 
increased scatter However, the likelihood fit is much less af- 
fected by large-scale structure than the directly computed me- 
dian redshift, which in contrast under-estimates the slope of the 
Zm(> ^606) r elation for Zm S 24.7 (see Fig.fTSll. This is the reason 
why IhoI who use the median redshift computed from spectro- 
scopic data in the magnitude range 21.8 < mm < 24 A, derive 
a significantly flatter Zm('M606) relation 
7H05 



-3.132 + 0.164OT606 (21.8 < meoe < 24.4) 



(43) 



leading to an estimate of Zm = 1.0 + 0.1 for their shear cata- 
logue. 

In order to verify the applicability of dTO for our fainter 
shear galaxies, we also plot Zm('«606) in Fig.fTslcomputed from 



photometric redshifts for the HUDF dCoe et al.ll2006l) . finding 
a very good agreement. 

Using the parameters (a,/3,a,b) we construct the redshift 
distribution for all GEMS and GOODS galaxies in our shear 
catalogue from (l42l) . The resulting redshift distribution has a 
median redshift Zm(GEMS + GOODS) = 1.46 ± 0.02(0.05), 
where the statistical errors stem from the uncertainty of a and 
b. Systematic uncertainties might arise from applying ( 1411 1 
for galaxies up to 1.5 magnitudes fainter than the magnitude 
range used to determine the fit. Additionally, the photometric 
redshift errors used in the maximum likelihood analysis do not 
tak e catastrophic outlier s or systematic biases into account, but 



Grazian et al. (l2006h for a comparison to the spectroscopic 



see 

subsample. Furthermore the impact of the large-scale structure 
on the source redshift distribution will be slightly different 
for the whole GEMS field compared to the GOODS region. 
We estimate the resulting systematic uncertainty as Az ^ 0.1, 
yielding Zm(GEMS + GOODS) = 1.46 + 0.02(0.05) ± 0.10. 
The constructed redshift distribution is well fit with 
a magnitude independent distribution ( l40b with 
(a,/3,zo) = (0.537,1.454,1.832). 

Given that we derive the redshift parametrisation from the 
matched GOODS-MUSIC galaxies in the magnitude range 
21.75 < m606 < 26.25, while a low level of redshift incomplete- 
ness already occurs for m606 <: 25.75 (see Fig.[T3]l. we repeat 
our analysis as a consistency check using only galaxies with 
21.75 < m606 < 25.75 yielding a very similar redshift distribu- 
tion with Zm = 1 -44. We thus conclude that the low level of in- 
completeness does not significantly affect our analysis. 

For the brighter galaxies in our shear catalogue with 
S/N > 5,m606 < 27.0, the constructed redshift distribution is 
expectedly shallower with Zm = 1.37 + 0.02(0.05) + 0.08 and 
can well be fit with a magnitude independent distribution 
(Unil with (a,/3,zo) = (0.529, 1.470, 1.717). Using our redshift 
parametrisation we also estimate the median redshift for the 
IH05 shear catalogue yielding Zm = 1 -25 ± 0.02 ± 0.08. Here we 
estimate slightly lower systematic eiTors due to the lesser ex- 
trapolation to fainter magnitudes. 

In Sect.|7] we will use our derived redshift distribu- 
tion to constrain cosmological parameters marginalising over 
the statistical plus systematic error in Zm- Furthermore we 
will use this redshift distribution when we compare cos- 
mic shear estimates for the GEMS and GOODS data with 
theoretical models. The theoretical cosmic shear predictions 
shown in this paper are calculated for a flat ACDM cos- 
mology according to the three year WMAP-only best-fitting 
values for (Qa, O.^, Qb, /1, n.s) = (0.76, 0.24, 0.042, 0.73, 0.95) 
( Spergel et al. 2006) for different power spectrum normalisa- 
tions cTg calcul ated using the non- Unear correction to the power 
spectrum froml Smith et al. ( 2003 ). 

At this stage we use the parallel data to test our pipeline 
and search for remaining systematics, while presenting a 
cosmological parameter estimation in a future paper based 
on a larger data set. Given the inhomogeneous depth and 
data quality of the parallel data, this cosmological parame- 
ter estimation will require a thoroughly estimated, field de- 
pendent redshift distribution. For the purpose of comparing 
the different estimators for shear and systematics to the ex- 
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pected shear signal in the current paper, we apply a simpli- 
fied global redshift distribution estimated from the F775W 
magnitudes in GOODS-MUSIC catalogue. Similarly to the 
F606W data we apply our likelihood analysis to all GOODS- 
MUSIC galaxies with 22.0 < nijjs < 26.0 yielding best fitting 
parameters (a,j8, a, b) = (0.723, 1 .402, 0.309, 0.395), for which 
Zm - 0.9395zo- The upper magnitude limit was chosen due to 
a similar turn-off point of Zm('M775) as in Fig. [15] indicating 
redshift incompleteness. To account for the different extinc- 
tion in the parallel fields and the CDFS (A™''^ = 0.017 mag), 
we apply an extinc tion correction based on the maps by 
Schlegelet aLl(ll998h . 

Using the extinction corrected magnitudes of all F775W 
galaxies in the parallel data shear catalogue, we con- 
struct a redshift distribution with z,,, - 1.34, which can 
be fit with a magnitude independent distribution with 
{a,p,ZQ) = (0.746,1.163,1.191). 

6. Cosmic shear estimates and tests for 
systematics 

In this section we compute different cosmic shear statistics and 
perform a number of diagnostic tests to check for the presence 
of remaining systematics. For the GEMS and GOODS data the 
plots in this section correspond to the larger galaxy set with 
S/N > 4 including the faint galaxies which are stronger af- 
fected by the PSF. 

6.1. Average galaxy ellipticity 

For data uncontaminated by systematics the average galaxy 
ellipticity is expected to be consistent with zero. Any signifi- 
cant deviation from zero indicates an average alignment of the 
galaxies relative to the pixel grid. We plot the average cor- 
rected but not rotated (see Sect. 15. Il l galaxy ellipticity (e™) 
for each field in Fig.|20] Whereas the global average is essen- 
tially consistent with zero for the GEMS and GOODS data 
({e'P = -0.0004 + 0.0011, <eif°) = 0.0012 + 0.0011), the av- 
erage ej''°-component is significantly negative for the parallel 
data«ef ) = -0.0084 + 0.0015, {e'^°) = 0.0020 ± 0.0015) cor- 
responding to an average orientation in the direction of the y- 
axis. 

6.1 .1 . Could it be residual PSF contamination? 

There are different effects which could in principle cause 
such an average alignment: For example one could spec- 
ulate that our PSF fitting technique fails for the parallel 
data or that our implementation of the KSB+ formalism 
under-estimates the PSF anisotropy correction, e.g. due 
to neglected higher-order moments. Yet, the average cor- 
rected galaxy ellipticity is consistent with zero for the 
GEMS and GOODS data, while the average uncorrected 
ellipticity is significantly non-zero for both datasets (par- 
allel: <ei) = -0.0102 + 0.0012, (ez) = 0.0028 + 0.0012; 
GEMS +GOODS : <e i ) = -0.0090 ± 0.0009, 

(62) = 0.0045 + 0.0009). Therefore this explanation be- 
comes quite implausible, particularly as the average number 



of stars usable to derive the fit is higher for the parallel data 
(Fig.|3. 

To further test whether imperfect PSF correction could be 
the cause, we plot the mean galaxy ellipticity as a function of 
the mean PSF anisotropy kernel on a field-by-field basis for 
parallel data in Fig. [16] While there is a substantial correlation 
between and the mean uncorrected ellipticity (ea) (corre- 
lation cor = cov[(^„), {ea)]/(o-(q^)ar(e„)) = 0.38), the mean PSF 
corrected ellipticity (e™) is basically uncorrected with 
(cor = 0.08), clearly indicating that imperfect PSF correction 
is not the culprit here. 

We also plot the mean corrected galaxy ellipticity (e™) 
computed in ^„-bins in Fig. [17] The absence of a correlation 
both for the GEMS+GOODS data and additionally <e™> in 
the parallel data again confirms the success of the PSF cor- 
rection. For the parallel data a moderate correlation is ob- 
served between (e™) and q[, which at first sight might be in- 
terpreted as an indication for imperfect PSF anisotropy correc- 
tion. However, it is important to keep in mind that qa is posi- 
tion dependent. Hence, if a different position dependent effect 
causes the non-zero (e™) it will also mimic a dependence on 
qa- From Fig. [5] we find for example that highly negative val- 
ues for qi appear mainly near medial y-positions close to the 
gap between the two chips. Thus, the apparent correlation be- 
tween (e™) and qi shown in Fig. [IT] could also be caused by 
a different effect which acts most strongly near the chip gap, 
such as CTE degradation (see Sect. |6.1.2] i or artefacts due to 
bad columns (see Sect. l6.1.3] l. In this sense the field-by-field 
comparison shown in Fig. [16] is a better test for imperfect PSF 
anisotropy correction, as it is independent of a possible posi- 
tion dependence. Given the fact that this test does not show a 
significant indication for imperfect PSF anisotropy correction, 
we conclude that it is most likely not the explanation for the 
non-zero (e™). We investigate the position dependence further 
in Sect. 16. 1.2] and compute the star-galaxy cross-correlation as 
an additional test for PSF anisotropy residuals in Sect. 16. 31 



6.1 .2. Impact of CTE degradation 

Another possible explanation is a degradation of the charge- 
transfer efficiency (CTE) due to charge traps created by the 
continuous cosmic ray bombardment. The ACS/WFC has two 
read-out amplifiers per chip, which are located in the four cor- 
ners of the instrument. Of major concern is the degradation 
of the parallel CTE causing charge trails behind objects in the 
readout direction, which also in the drizzled images is approx- 
imately parallel to the y-direction. These charge trails lead to 
an average alignment of objects in the y-direction, correspond- 
ing to a negative average ei ellipticity component. As the depth 
of charge traps is limited, faint objects loose a larger fraction 
of their charges than bright ones, leading to a strong signal-to- 
noise dependence of the effect. Therefore the PSF correction 
estimated from high signal-to-noise stars does not provide a 
sufficient CTE correction for faint galaxies. 

For a uniform distribution of charge traps the impact of 
CTE degradation depends linearly on the number of parallel 
transfers, so that objects located near the gap between the two 
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Fig. 16. Mean galaxy ellipticity before 
«e„), left) and after «e™), right) PSF cor- 
rection as a function of the mean PSF 
anisotropy kernel averaged over all galax- 
ies in a field {qa), computed on a. field-by- 
field basis for the F775W parallel fields. 
The lack of a correlation after PSF cor- 
rection (cor = 0.08) is a clear indication 
that PSF anisotropy residuals cannot be 
the origin for the negative average ellip- 
ticity <ej'°). 
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Fig. 17. Mean PSF corrected galaxy el- 
lipticity <ej^°> binned as a function of the 
PSF anisotropy kernel for the paral- 
lel data (left) and the GEMS+GOODS 
data (right). The binning (indicated by 
the horizontal error-bars) was chosen such 
that all bins contain an equal number of 
galaxies. The lack of a correlation for the 
GEMS+GOODS data and {e'^") for the 
parallel data confirms the success of the 
PSF correction. The interpretation of the 
moderate correlation detected for (e™) in 
the parallel data is ambiguous as it can 
also be caused by a position dependence 
of<ef). 
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Fig. 18. Average corrected galaxy ellipticity {e™) for the paral- 
lel F775 W galaxy fields as a function of Ay, the y-position rela- 
tive to the gap between the two camera chips. The curve shows 
{e^^°){h.y) box-averaged over 3000 galaxies. For certain Ay the 
error-bars indicate the width of the averaging in Ay and the 
error of the estimate. The straight lines indicate the expected 
dependence if the negative (e™) was purely caused by CTE 
degradation assuming a linear dependence of the mean elliptic- 
ity on the CTE charge loss. 
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Fig. 19. Average coiTected galaxy ellipticity (e™) for the paral- 
lel F775W galaxy fields as a function of time since the installa- 
tion of ACS on March 7, 2002. The solid line shows a linear fit. 
If the negative (e™) would be created by degradation of CTE 
an increase of the effect with time would be expected, which is 
not supported by the data. 




chips will be affected the most. lMutchler & Siriannil (120051) find 
no significant difference in the parallel CTE for the two chips. 
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indicating that also any impact on the weak lensing measure- 
ment should be symmetric between the two chips. In Fig.fTSl 
we plot (e''^°) as a function of Ay, the y-position relative to the 
gap between the two camera chips. Although for the lower chip 
(Ay < 0) (e'^°)(Ay) roughly agrees with the linear trend ex- 
pected for a CTE degradation, there are significant deviations 
for the upper chip (Ay > 0). 

Furthermore the ACSAVFC C TE decr eases nearly 



linearly with time (iRiess & Mackl 120041; iRiessI 12004 



Mutchler & Srrianni l2005h so that one would also expect 
a Unear decrease of (e™) with time, which is not in agreement 
with the data (Fig.[T9Tl. In addition, again, the discrepancy does 
not occur for the GEMS and GOODS data, which were taken 
nearly in the same time period as the parallel data. We thus 
conclude that CTE degradation is not the domi nant source for 
the o bserved negative (e™). Note that Rh odes et al.l (12005 



20071) detect discrepancies between their focus-dependent 
TINYTIM model and stars in the COSMOS data, which they 
interpret to likely be caused by CTE degradation. This is not in 
contradiction to our results, as the COSMOS data were taken 
at later epochs with significantly increased CTE degradation. 

As a further test we also split the data shown in 
FiguresfTSl and [19] into a low and a high signal-to- 
noise subset. Here the observed dependencies are qual- 
itatively unchanged, but at a lower significance, with 
a slightly larger absolute values of the negative {e^^°} 
for the fainter sample: <e\'°>(S/N < 7.5) = -0.0092 ± 0.0022, 
{e'p(S/N > 7.5) = -0.0077 + 0.0018. If the effect was caused 
by CTE, one would probably expect a stronger dependence on 
the signal-to-noise ratio. 



6.1.3. Impact of dithering 

In order to understand the origin of the negative (e™) for the 
F775W parallel data it is helpful to consider the differences be- 
tween the two surveys, as the problem does not occur for the 
F606W GEMS and GOODS images. Besides the different fil- 
ters and more homogeneous depth of the GEMS and GOODS 
tiles there are only two effects which can significantly affect 
the image quality: Firstly the F775W fields are taken in par- 
allel in contrast to the F606W data. Although this could have 
some impact on the image PSF (Sect. l3.ll l. it is taken into ac- 
count in our PSF correction scheme (Sect. 14.4] ). Secondly the 
GEMS and GOODS data are well dithered, whereas most of 
the parallel fields were observed with no or only small dither- 
ing as defined by the primary observations. To test the impact 
on the galaxy shape measurement we split the parallel fields in 
Fig. |20] into three sets according to the maximal shift between 
the exposures in jc-direction AX. Indeed (e™) is almost consis- 
tent with zero for the well dithered fields with AX > 10 pixels 
(6 WFC pixels), whereas it is significantly negative for the less 
dithered fields. 

Proper dithering is important to correct for bad or hot pix- 
els, which otherwise create artifacts in the co-added frame. 
Without dithering known bad pixels lead to output pixels re- 
ceiving zero weight, which we set to zero pixel value, while un- 



known bad pixels such as spontaneously hot pixels or variable 
bias structures directly contribute with their bad pixel value. 

Bad pixels are not completely randomly distributed on the 
CCD chips, but sometimes occur as bad columns or clusters 
of bad pixels, which are preferentially aligned in the readout 
direction and therefore the y-direction. Thus, without proper 
dithering the shapes of faint objects containing bad columns or 
pixel clusters could possibly be influenced such that a slight 
average alignment in the y-direction is created and a negative 
(e'j^°) is measured. We expect that faint galaxies are stronger, 
and due to their size more likely, affected than compact high 
signal-to-noise stars, which additionally might be rejected as 
noisy outliers during the PSF fitting, explaining why this effect 
is not taken into account by the PSF correction. 

We try to minimise the impact of known bad pixels by 
rejecting galaxies containing low weight pixels within their 
SExtractor isophotal area (see Sect. l3.3l l. However, also a bad 
column located near the edge of a galaxy image might bias the 
shape estimate without being rejected in this way. Note that bad 
column segments appear with a higher density near the chip 
gap, which might qualitatively explain the Ay dependence plot- 
ted in Fig.[Il 

Although the comparison shown in Fig. |20] supports our in- 
terpretation that the negative (e™) is caused by a lack of dither- 
ing, we will need to further investigate this effect on the basis 
of the complete ACS Parallel Survey for a final judgement, as 
it extends over a much larger time span allowing a clearer dis- 
tinction from CTE effects. 

So far we co-add parallel data observed within one visit 
to maximise the stability of the image conditions. Due to the 
successful PSF correction for the two-epoch GOODS data (see 
also Sections [6.3! and l6r4] i we are confident that a combination 
of different visits will also be possible for parallel data with re- 
observations, which will reduce the number of fields with poor 
dithering. Additionally we are working on an improved search 
algorithm for galaxies which are affected by bad columns. 



6.2. Error estimates 

In the following subsections we compute several estimators for 
the cosmic shear signal and remaining systematics. The statis- 
tical errors of these estimates are always computed in a similar 
way. 

6.2.1. Parallel data 

Bootstrapping on galaxy basis. To derive statistical weights 
for ^± and (M^p), we generate for each field / 200 bootstrap 
samples of the galaxy catalogue and compute ^±jj and {M^^)ij 
for each angular bin j. The weight w,j for this field and bin is 



then given as the inverse bootstrapping variance w,y 
yielding the combined estimates 



fields J. 



yi ^fields 



I/CT 



(44) 



Wij 
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Fig. 20. Average corrected galaxy elliptic- 
ity for the parallel data (left) and the 
GEMS+GOODS data (right). The open 
symbols represent single field averages, 
whereas the bold symbols with error-bars 
(Icr) show global averages. The large 
error-bars in the lower right comer cor- 
respond to the average single field Icr 
error, where the GEMS+GOODS error- 
bar is smaller compared to the paral- 
lel data error-bar due to the higher av- 
erage galaxy number density (Fig.lllt. 
The parallel data was split according to 
the maximal dither between the expo- 
sures in the x-direction AX [pixels], as 
the x-dithering determines the possibili- 
ties to correct for bad columns. 



The estimate for the galaxy-star cross-correlation (see 
Sect. 16. 31 ) is calculated accordingly, with bootstrapping of the 
galaxy catalogue and a fixed stellar catalogue. 

Bootstrapping on field basis. We determine the measurement 
error of the field combined estimates for ^j-j and {Ml^)j from 
300 bootstrap samples of our fields, combining the estimates 
for each realisation according to ( l44l i. The error of the com- 
bined signal in each angular bin j is then given by the bootstrap 
variance cr^. This error estimate accounts both for the shape 
noise and cosmic variance. 



the jackknife estimate for is given by the average 
l-t j - (^* .j), and the jackknife estimate of the covariance be- 
tween bins i and k can be computed as 

<A^.,A^,,) = j^^^ g (^:,,, - 1,,,) (^:,,, - 1,,,) . (46) 

Note that this jackknife method is expected to slightly under- 
estimate the error even on scales much smaller than the field 
size due to the mixing of power between different scales in the 
non-linear regime. 



6.2.2. GEMS and GOODS 

Bootstrapping on galaxy basis. For the combined GEMS and 
GOODS mosaic catalogue we analogously perform bootstrap- 
ping on galaxy basis to derive the shape noise error. The eiTors 
plotted for the galaxy-star cross-correlation and the E-/B-mode 
decomposition within Sections |63] and 16.41 correspond to this 
bootstrap variance. 

For the cosmological parameter estimation in Sect.|7] co- 
variances are required, which additionally take sampling vari- 
ance into account. We compare covariances estimated directly 
from the data using a jackknife method with estimates from 
Gaussian realisations of the cosmic shear field. 

Jackknife method. We use the modified jackknife method ap- 
plied by H05 to estimate the covariance matrix of the cosmic 
shear estimators. In contrast to the bootstrapping on galaxy 
basis, the jackknife method applied includes an estimate for 
small-scale cosmic variance. However, it must under-estimate 
cosmic variance on scales of the order of and larger than the 
field size. We describe the algorithm in terms of the correlation 
functions ^j-: We first compute the correlation function ^j- y in 
the angular bin j from the complete galaxy catalogue. Next, 
we divide the whole survey into separate sub-regions on the 
sky, where for convenience we use the A^ = 78 individual ACS 
tiles. Then, the correlation function ^-t ,j is computed omitting 
the /-th subregion for i - 1, A. With 

C7 = ^^±.i-(^-l)^±.'V' (45) 



Sampling variance from Gaussian random fields. Given that 
the GEMS and GOODS mosaic samples only one particular 
field in the sky, the large scale sampling variance errors cannot 
be determined from the data itself. In order to derive a theo- 
retical error estimate we have created 2000 1° x 1° Gaussian 
realisations of the shear field for a ACDM cosmology with 
cTg - 0.7 and the GEMS redshift distribution, which we pop- 
ulate with 96 galaxies arcmin"^ with ellipticities randomly 
drawn from our shear catalogue. We then select a ~ 28' x 28' 
subregion representing the actual masked geometry of the mo- 
saic. From the sheared ellipticities we then compute the co- 
variance matrix of the correlation functions from the different 
realisations (see lSimon et al ]|2004 . This provides us with a ro- 
bust estimate of the error covariance in the Gaussian limit also 
including the shape and shot noise contribution. Note, how- 
ever, that the Gaussian assumptio n strongly under-estimates the 
sampling varianc e for e < 10' dKilbinger & Schneideij|200l 
Semboloni et al ] |2007h . which we further discuss in Sect.|2l 



6.3. Star-galaxy cross-correlation 

An important diagnostic test for the effectiveness of the PSF 
anisotropy coiTection is given by the cross-correlation between 
uncorrected stellar ellipticities e* and coiTected galaxy ellip- 
ticities y, which can be used as an estimate for residual PSF 



contamination. Following Ba con et al.l (120031) we compute 



{ye*)m{7e*)m 
{e*e*}(0) 



(47) 
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For the parallel data we substitute e* with the smearing cor- 
rected PSF model ellipticity 



^mod,Qr 



^ 2Ccal sm DRZ (-^ ^ „ ^ 
gal 



(48) 



(see Eg. [39]), at all galaxy positions, which is necessary as 
is very noisy and undetermined in many bins due to the 
few stars present in most of the single parallel fields. 

As can be seen from Fig. [21] C^^^ is consistent with zero for 
the GEMS and GOODS data for all B indicating that the PSF 
correction works very well for this dataset. For comparison we 
also plot C^^^ computed from the smearing but not anisotropy 
corrected galaxy ellipticities, which exceeds the theoretically 
expected cosmic shear signal, emphasising the need for proper 
PSF correction. 

In contrast, C^'^^ is non-zero for the parallel data for most 
Q. Considering the results from Sect. l6.1.3l we interpret this re- 
maining systematic signal as cross-correlation between the (av- 
erage) PSF pattern and the mean ellipticity component induced 
by the lack of dithering. This interpretation is supported by the 
fact that C^^^ is almost consistent with zero when computed 
from the corrected galaxy ellipticities minus the mean elliptic- 
ity (Fig.lSTTl, suggesting that the PSF correction also performs 
well for the parallel data. 

The underlying assumption of our PSF correction algo- 
rithm is that the stellar fields sample the parameter space of PSF 
variations in the galaxy fields sufficiently well (see Sect. |4.4.3] l. 
To test this assumption we repeat the analysis always using the 
iecont/-best fit PSF model instead of the best fitting model. 
If the sampling of the PSF variations was not sufficient, we 
would expect a significant impact on the PSF corrected el- 
lipticities and particularly C^^^ when switching to the sec- 
ond-hest fit PSF model. However, as the observed impact is 
negligible both for C^^^ (left panel of Fig.l2TI) and the mean 
corrected galaxy ellipticity ({e;'"'™''^) = -0.0085 ± 0.0014, 
^^i^cmodiy ^ Q QQ^g ^ 0.0014, compare to Sect.|6l]l, the sam- 
pling of the PSF parameter space indeed seems to suffice. 

6.4. E-ZB-mode decomposition 

As a further test for contamination of the data with systemat- 
ics we decompose the shear signal into E- and B-modes us- 
ing the shear correlation functions ^^(0), ^^(0) (Fig.l22l) and 
the aperture mass dispersion (Fig. [23]). For this we first calcu- 
late ^+(0) and ^_(0) in 300 (1800) finite linear bins of width 
A6' = a.'83 (1'.'17) from l"to 4.'2 (35') for the parallel (GEMS 
and GOODS) data, ^^-^(fl) and {MI^ J{0) are then computed 
according to equations ( 18115116b and logarithmically re-binned 
to reduce noise. 



6.4.1 . ^^/^^ decomposition 

As the computation of ^^'^(0) requires knowledge of ^_ also 
for 6 larger than the field size (see Sect. l2.1l ), we substitute the 
measured ^_ for 6 > 4' (9 > 35') with theoretical predictions 
for a fiducial ACDM cosmology with erg = 0.7. The impact of 
the fiducial cosmology on the E-/B-mode decomposition can 



be estimated by comparing ^^'^(9) computed for erg = 0.6 and 
CTg = 1.0 (dotted lines in Fig. 122b. While the difference is small 
for the GEMS and GOODS data (~ 2 x 10"^), the small size 
of the single ACS fields leads to a stronger cosmology depen- 
dence (~ 1.5 X 10 "*) for the parallel data. The B-mode compo- 
nent is consistent with zero for both datasets indicating that 
we are not subject to major contaminations with systematics. 
The only exception is the slightly negative for the GEMS 
and GOODS data at large scales, which is an artefact of the 
discontinuity between the fiducial cosmological model and the 
low shear signal measured at large scales (see the E-mode sig- 
nal and Sect. l6.5.T]) in combination with the bootstrap errors, 
which do not take cosmic variance into account. 

6.4.2. {mI^)/{MI) decomposition 

Also the B-mode component of the aperture mass dispersion 
{M]){9) is consistent with zero for both datasets indicating 
the success of our PSF correction scheme (Fig.l23ll. Note that 
the E-/B-mode mixing due to incomplete kno wledge of 4^+ (6*) 



for small 0, which was recently discussed by iKilbinger et al 
(2006), only leads to minor effects for the 9 range consid- 
ered here, since we truncate ( +{0) only for < 9^^^ - 2". See 
Schneider & Kilbingeii ( 120071) for a E-/B-mode decomposition 
which can also be used for larger flniin- 

6.5. Shear correlation functions 

We plot our estimate for the logarithmically binned shear two- 
point correlation functions {yiyt){0) and {yx7x)(0) in Fig.l24l 
Note that we use ^±(0) for the cosmological parameter estima- 
tion in Sect.|7] but plot the equivalent data vector s (y>y >)(0) and 



(TxTx)(S) in order to enable the comparison with|H05 



6.5.1 . GEIVIS and GOODS data 

As we have shown in the previous sections, the GEMS and 
GOODS data are not contaminated with significant non-lensing 
signals. We are therefore confident that the measured shear sig- 
nal (right panel of Fig.l24]i is of cosmological origin. While 
we detect significant shear correlations at small angular scales 
consistent with predictions for erg ~ 0.6, both {ytyt){0) and 
{yxyx){0) are consistent with zero for 9 > 5', which we inter- 
pret as caused by a large-scale under-density of the foreground 
structures in the CDFS. 

There is good agreement between the error-bars determined 
from the jackknife method and from Gaussian realisations. 
Only for scales of the order of the field size the jackknife 
method significantly under-estimates the modelled errors as it 
does not account for large-scale cosmic variance. Note the good 
agreement of the data with the results from lH05i 

6.5.2. Parallel data 

While the measured shear correlation functions are roughly 
consistent with the plotted ACDM predictions for erg ~ 0.8 
(left panel of Fig.l24li. one must be careful with its interpre- 
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Fig. 21. Star-galaxy cross-correlation for the parallel data (left) and the GEMS and GOODS data (right), where C'>" is 
calculated from the uncorrected stellar ellipticities e* for the GEMS and GOODS data and the PSF anisotropy model e*^^ for the 
parallel data. The squares show C'^^'^ computed from the corrected galaxy ellipticities. For the parallel data this can be compared 
to the crosses (stars), where the PSF correction was derived using the second-best fit PSF models. The negligible difference 
between the two indicates that the F775W stellar field exposures sample the PSF variations sufficiently well. For comparison we 
also plot C^^^ determined from the smearing but not anisotropy corrected galaxy ellipticities (triangles), and in case of the parallel 
data also computed from corrected galaxy ellipticities after subtraction of the mean corrected ellipticity (circles). The different 
data sets are displayed with different 0-offsets for clarity. The dashed (dotted) line shows ACDM predictions for {jiy,} ((yxTx)) 
for erg = 0.7. 



tation due to the detected indications for remaining systematics 
(Sect. l6.l] and 16.31) . even if they do not show up as B-modes. 
We thus postpone the cosmological interpretation of the paral- 
lel data shear signal to a future paper based on a larger data set 
with further corrections for the remaining systematics. 

7. Cosmological parameter estimation from the 
GEMS and GOODS data 

Having shown that our GEMS and GOODS shear catalogues 
are not subject to significant non-lensing systematics, we use 
our estimate of the shear correlation functions, binned in 14 
logarithmic bins for 0.'058 < 6 < 28.' 1, in combination with the 
determined redshift distribution (Sect. 15. 3] ) for a cosmologi- 
cal parameter estimation usin g a Monte Carlo M arkov Chain 



( MCMC) technique ( s ee e.g . iTereno et al.l 120051) as detailed 
in iHetterscheidt et alj ( 120061) . Here we utiUse the covariance 



matrix derived from the Gaussian realisations. This is moti- 
vated by the good agreement with the errors determined from 
the jackknife method at small scales indicating rather low 
impact of non-Gaussianity. However, using ray-tracing sim- 
ulation s iKilbinger & Schneideii (l2005h and ISemboloni et al. 
(l2007h found that Gaussian statistics strongly under-estimate 
the covariances also for GEMS like surveys, which we further 
discuss below. 

For the parameter estimation we consider two simple 
ACDM cosmological models: 

A: a A-universe with Q^, Qa 6 [0, 1.5], 



B: a flat universe: Qm + = 1 with > 0, 

both with fixed {w,Q.h,ns) - (-1,0.042,0.95). We assume a 
strong constraint h = 0.70 ± 0.07 for the Hubble param- 
eter, as supported by the HST key project (jpreedman et al] 
200 Ih and co mpute the non-line ar power spectrum us- 



ing halofit (iSmith et alJ l2003j). with t h e sha pe pa- 



rameter calculated according to [sugiyamal ( 1995 ): F - 
Qm/i ex p [-Qb(l + V2/i/^m)L and the transfer function as 



given in Efstathiou et al.l ( 11992 ). In the likelihood analysis we 



marginalise over the uncertainty in both h and our redshift dis- 
tribution. 

We plot the derived likelihood contours for 0-^,0.^, 
and Qa in Fig.|25] where we use all galaxies with S/N > 4 
(A^gai = 96 arcmin"^) corresponding to a median redshift 
Zm = 1.46 ± 0.12. For the more general model A the data only 
weakly constrain QA(Qm - 0.3) - 0.64^q^j, whereas more 
stringent constraints are found for cr8(Qni - 0.3) - 0.52^q [^, 
or respectively, Qm(o'8 = 0.6) - 0.26^QQg, reflecting the 
marginalised 68% confidence regions with strong priors on 
Q,n or cTg respectively. Assuming flatness (model B) changes 
the estimates only marginally to cr^i^m = 0.3) = 0.5 1^" ?? and 



n^io-s = 0.6) = 0.25!j;°^, respectively. 

Using the more conservative sample 
S/N > 5, OT606 < 27.0, A^gai = 72 arcmin"^, 
leads to a higher estimate of crg(Q 



-0.13 

selection with 
Z,n = 1.37 ±0.10 
0.3) = 0.59+°" 



-0.14 



or 



Q,„(cr8 - 0.6) - 0.30^qq|| for model A without significantly af- 
fecting the error. In principle, one would expect that the inclu- 
sion of the faint galaxies increases the signal-to-noise of the 
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Fig. 22. E-/B-mode decomposition of the correlation functions for the parallel data (left) and the combined GEMS and GOODS 
data (right). The open circles show and computed using a fiducial ACDM model with erg - 0.7 for the extrapolation in 
Eq. (|9]l, whereas the dotted lines correspond to erg, - 1. (upper line for f^, lower line for and erg - 0.6 (lower line for^f, 
upper line for ^^). The thin solid lines show ACDM predictions for erg = (0. 6, 0.8, 1.0). In the right panels we also plot the IH05 
GEMS only estimate for and for erg = 0.7 (crosses). Note that the |H05| catalogue is slightly shallower 
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Fig. 23. E-/B-mode decomposi- 
tion of the aperture mass disper- 
sion for the parallel data (left) 
and the combined GEMS and 
GOODS data (right). The thin 
solid lines show ACDM predic- 
tions for erg = (0.6,0.8,1.0). 



shear measurement as both the galaxy number density and the 
lensing efficiency increase. However, we can confirm the trend 
seen bv IhOSI that the faintest galaxies appear to mainly add 
noise and dilute the signal. This is also consistent with the re- 
sults from the STEP2 image simulations, where we find that 
the shear calibration of our KSB+ implementation is on aver- 
age accurate to ~ 3%, but shows a significant dependence on 
magnitude, with a slight over-estimation at the bright end and 
a ~ 20% under-estim ation of the shear for the faintest galaxies 
dMassev et al. I l2007ah . Given the on average good calibration 
found for our analysis of the STEP2 simulations, which incor- 
porate a cut S/N > 4, we consider the estimate of erg for the 
same cut to be more robust. Yet, as the magnitude and size dis- 
tribution, and additionally also the noise correlations are some- 
what different for the STEP2 simulations and our data, we ex- 



pect a slight remaining systematic error also for the average 
shear calibration. Therefore, we use the difference of the two 
estimates for erg as a conservative estimate of this uncertainty 
and take it into account as additional systematic eiTor, yielding 
o-g(n„ = 0.3) = 0.52+JJ [i(stat) + 0.07(sys). For a future shear 
tomography analysis this issue will need to be revisited, as it 
does not only require accurate shear calibration on average, but 
also over the whole magnitude range. 

The constrained value for erg is significantly lower 
than the estimates from other recent le nsing sur- 



veys, e.g. 
erg = 0.85 : 



erg = 

0.06 



86 ± 0.05 dSemboloni et alj l2006l) 



dHoekstra et al 



2006h. 



both from 



and 
the 



CFHTLS for Q.^ = 0.3. see i Hetterscheidt et"ai] (12006') for a 
compilation of recent estimates. Our results are consistent with 
erg = 0.8, Qm = 0.3 only at the 3cr-level assuming Gaussian 
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Fig. 24. Two-point correlation functions (7,7,) and (yxTx) for the parallel data (left) and the combined GEMS and GOODS data 
(right). In the right panels we plot our estimate (open circles) both with the er rors d etermined from Gaussian realisations (bold 
error-bars) and the Jackknife errors (thin caps), and for comparison also the IHOSI results (crosses). The thin solid lines show 
ACDM predictions for erg - (0.6,0.8, 1.0). Note the very low cosmic shear signal measured from the GEMS and GOODS data 
for large 9. 



cosmic variance, which we interpret as a substantial local 
under-density of the foreground structures in the CDFS. 
In order to allow a clear comparison to the H05 results, 
who determine cr8(Qn,/0.3)*^'*^ =0.68 + 0.13, we recompute 
our redshift distribution using their Zin(mag) relation (|43] |. 
yielding a median redshift Zm - 1-12 (Zm = 1.07) for the 
galaxies with S/N>4 (S/N > 5, msoe < 27.0). Then we 
repeat the cosmological parameter estimation as sumin g a 
redshift uncertainty Azm = 0.1 to be consistent with HOS". For 
this redshift distribution we find cr^in^ - 0.3) - 0.62+° 
(cr8(Q m = 0.3) - 0.66^q J4) in excellent agreement with the 
5I 



IHO5 I results. We thus conclude that our lower erg estimate 
compared to H05 is mainly a result of our new redshift 
distribution based on the GOODS-MUSIC sample, and that 
the two independent shear pipelines yield consistent results 
(see also Sect. |5.2t . 

Our estimate of the statistical eiTor includes the shape 
noise contribution, the estimated uncertainty of the redshift 
distribution, and a Gaussian estimate for cosmic variance. 
Although there is good agreement of the eiTors from the jack- 
knife method and Gaussian realisations at small scales, we 
expect to under-estimate cosmic varian ce d ue to non-linear 
evolut ion. Kilbin ger & Schneider (I2OO5') and 'Semboloni et aT 



(I2007h found that the Gaussian approximation can lead to a 
substantial under-estimation for the coiTelation function co- 
variance matrix in the non-linear regime. Using a fitting for- 
mula found bv ISembo loni et al. (2007) we estimate that the 
diagonal elements of the ^+ covariance matrix will be under- 
predicted by a factor of ~ 2.9 for a single source redshift plane 
at z = 1.4. As this corresponds to the median redshift of our 
galaxies and since our shear signal has the highest significance 



forO'6 < 6 < 5' (see Fig.l24li. which (logarithmically averaged) 
roughly corresponds to a scale ~ 2', we estimate very broadly 
that we on average under-estimate the cosmic variance contri- 
bution to the covariance matrix by a factor of ~ 2.9 leading to 
an error of erg which is actually larger by ~ y/lS ~ 1 .7. 

Apart from the shear calibration uncertainty considered 
above, further systemati c errors might be introduced by intrin- 
sic ali gnment of sources dBrown et al .I2OO2H iKing & Schneiden 



'2002' Hevmans & Heave nsl2003l;lHevmans et al.ll2004l2006cl: 
Mandelbaum et al. 2006t) or a correlation between the intrin- 
sic elUpticities of galaxies and the densit y field responsible 



for gravitational lensing shear, detected by ' Mandelbaum et al 
(2006). Given the depth of the data analysed here, we how- 
ever expect that the impact of these two eff'ec ts will be small 



compared to the statistical uncertainties (see Hevm ans et al 



l2006c ). Further uncertainties arise from the limited accuracy of 
the predictions for the non-linear power spectrum. Yet, given 
that the measured shear signal is particularly low for large 
(see Sect. l6.5.l1 l. which are less affected by non-linear evolu- 
tion, this cannot explain the low estimate of erg for the GEMS 
and GOODS data. 

8. Conclusions and outlook 

We have presented a cosmic shear analysis of a first set of 
HST/ACS pure parallel observations and the combined GEMS 
and GOODS data of the CDFS. We estimate that our new cor- 
rection scheme for the temporally variable ACS PSF reduces 
the systematic contribution to the shear correlation functions 
due to PSF distortions to < 2 x 10"^ for galaxy fields contain- 
ing at least 10 stars. This is currently the only technique tak- 
ing the full time variation of the PSF between individual ACS 
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Fig. 25. Constraints on erg, Qj^, and Qa from the GEMS and GOODS data using all galaxies with S/N > 4. The three large contour 
plots show marginalised joint 2-dimensional 1,2, and 3cr likelihood contours for model A, whereas the small one was computed 
assuming flatness (model B). For model A the marginalised probability is plotted on the right for Q.^ (top), Qa (middle), and erg 
(bottom), where the thick solid curves correspond to the total marginalised values, while the thin solid (dotted) Unes correspond, 
from top to bottom, to fixed Qa = 0.7 (erg = 0.6), Q™ = 0.3 (erg = 0.6), Qa = 0.7 (Q^ = 0.3). 



exposures into account. In the GEMS and GOODS data the 
success of the PSF correction is confirmed by a number of 
diagnostic tests indicating that the remaining level of system- 
atics is consistent with zero. For the parallel data we detect 
a low level of remaining systematics manifesting in a slight 
average alignment of the measured galaxy ellipticities in the 
y-direction, which we interpret to be due to a lack of proper 
dithering. We are currently further investigating this effect and 
exploring ways to correct for it, which will be necessary for 
the cosmic shear analysis of the complete set of ACS paral- 
lel observations. Although the degradation of the ACS charge- 
transfer-efficiency has not been found to be a problem for the 
early data analysed in this work, an in-depth analysis and cor- 
recti on will probab ly be required for the complete data set (see 
also [Rhodes et al.li2005, 2007) . Furthermore the parallel data 
are rather inhomogeneous regarding depth and extinction, rais- 
ing the need for a well calibrated field-dependent redshift dis- 
tribution. It will also be necessary to carefully exclude any se- 
lection bias which might arise for certain classes of primary tar- 
gets, particularly galaxy clusters. Once these remaining obsta- 
cles are overcome, it will be possible to measure cosmic shear 
at small angular scales with unprecedented accuracy from the 
complete ACS Parallel Survey, with a strong reduction both 
of the shape noise and cosmic variance error due to many in- 



dependent pointings. The main limitation of the cosmological 
interpretation of the data might then arise from the current ac- 
curacy of theoretical predictions for the non-linear power spec- 
trum at small scales. An interestin g comparison will be possi- 
ble with the ACS COSMOS data dMassev et ai]|2007 bV from 



which cosmic shear can be measured on a wide range of angu- 
lar scales. 

Given the high demands concerning the control over 
systematics for cosmic shear measurement s with ACS, the 
derived technical expertise (see also HOSt Jee et al. _ 2005 
Rhodes et al1l2005l: iRhodes et al.ll2007t iLeaufliaud et al.ll200 



will also be of benefit for other weak lensing studies with the 
instrument, and possibly also other research fields requiring ac- 
curate PSF modelling. 

Due to the weakness of the shear signal on the one hand, 
and the strong impact of poorly understood systematics on 
the other hand, an analysis of identical datasets with more 
than one independent pipeline is of great value to check 
the reliability of the algorithms employed. In this work we 
have independently re-analysed the ACS observations of the 
GEMS and GOODS fields. If we assume the same redshift 
parametrisation, our shear estimates are in excellent agreement 
with the earlier results found by H05 indicating the reliabil- 
ity of both lensing pipelines. Such an independent compari- 
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son will also be highly desired both for the complete ACS 
Paral lel Survey (Rh odes et al. in prep.) and the ACS COSMOS 
field (IMassev et al.i,2007b.) . These comparisons, together with 
the results from the STEP project, will aid the preparations 
of future space-based cosmic shear survey such as DUNE 
or SNAP, which will reach a very high statistical accuracy 



dRefregier et al.1 120041) requiring the conti nued advancement 
of irn p roved algorithms such as shapelets ('Bernstein & Jarvis' 
Refreder & BaconI 120031: Massev & Refregier. .2005; 



M . -| | - — if-T' 

Kuijkenll2006l:lNakajima & Bernsteinll2007h . 



Finally, we want to stress the possible impact of the field 
selection on a cosmic shear analysis: The Chandra Deep Field 
South was originally selected in a patch of the sky charac- 
terised by a low Galactic neutral hydrogen c olumn density 



(A^H = 8 X 10 cm ) and a lack of bright stars dGiacconi et al 



20011) . Additionally, it neither contained known relevant extra- 



galactic foreground sources nor X-ray sources from the ROS AT 
ALL-Sky Survey Catalogu ^ excluding e . g. the presence of a 



low redshift galaxy cluster. lAdami et al.l (120051) present a de 



tailed analysis of compact structures in the CDFS showing 
the presence of a chain-like structure at z = 0.66, a mas- 
sive group at z = 0.735 embedded into a galaxy wall ex- 
tending beyond t he 21' x 21' field covered by the Vimos 
VLT Deep Survey (ILe Fevre et al.,,2004). and a further massive 
group at z = 1 .098 (see a lso lOilli et al 



20031: ISzokolv et al 



2004i:IVanzella et al.ll2006l) . IWolf et al.1 i2004) identify a strong 



galaxy over-density at z ~ 0.15, which is too close to produce 
a significant lensing signal. Given the lack of massive struc- 
tures at lower redshifts 0.3 < z < 0.6 with high lensing ef- 
ficiency, one would expect to measure a shear signal biased to 
lower values in this field as a result of strong sampling variance. 
Therefore it is not surprising that our local single field esti- 
mate of crs,cDFs(nm = 0.3) = 0.52+°;[ '(stat) + 0.07(sys) based 
on a source reds hift distribution der ived from the GOODS- 
MUSIC sample (Grazian et al.' '2006^, is incompatible at the 
3cr-level assuming Gaussian cosmic v ariance with recent re- 



sults of other weak len sing studies (e.g. iHoekstra et al. I l2006t 



ISem boloni et al.' '2006"), which probe much l arger regions on 
the sk y. Kilbinger & Schneider (20051) and ISemboloni et al 



(l2007h investigate the impact of non-Gaussianity on cosmic 
shear covariances. From their results we broadly determine an 
under-estimation of the cosmic variance contribution to our er- 
ror on cTs by a factor s; 1 .7, indicating that the CDFS is still 
an exceptionally under-dense field, but with a lower signif- 
icance (~ 2cr) than under the Gaussian assu mptio n. Our erg 
estimate is also significantly lower than the IH05| results of 
cr8(Q„/0.3)°-^^ = 0.68 ± 0.13 due to the deeper redshift dis- 
tribution found in our analysis with a med i an so urce red- 
shift Zm = 1.46 ± 0.12. Recently IPhleps et all (l2006h found a 
strong deficiency of faint red galaxies in the CDFS for the 
redshift range 0.25 < z < 0.4 indicating a substantial under- 
density, which is in excellent agreement with the low shear sig- 
nal found in our analysis. 

We believe that the CDFS represents a somewhat extreme 
case. However, also other cosmic shear studies which observe 
a low number of small "empty fields" could be slightly bi- 



ased just due to this prior selection. Such a bias can of course 
be eliminated either with the observation of sufficiently large 
fields or truly random pointings, which are realized in good ap- 
proximation for a large fraction of the fields in the ACS Parallel 
Survey. 

We plan to further investigate the peculiarity of the CDFS 
based on a shear tomography analysis with photometric red- 
shifts derived for the full GEMS field, also using deep ground- 
based _02tical_jmages_jTom the MPG/ESO 2.2m telescope 



dHildebrandt et al. I l2006h m combination with i nfrared images 



from the ESO 3.5m NTT dOlsen et al.ll2006blah . If the low es- 
timate for cTg cQps indeed stems from an under-density of fore- 
ground structures we would expect an increased shear signal 
for a high-redshift sample of source galaxies due to the spec- 
troscopically confirmed structures at z = 0.735 and z = 1.098. 
Comparing the results with ray-tracing through N-body simu- 
lations we aim to further quantify the rarity of such an under- 
dense foreground field. 
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